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ABSTRACT 


A  theoretical  and  experimental  study  of  signal  propagation  through  a  series  of 
cascaded  rectangular  cavities  connected  by  walls  containing  narrow  slots  and  with  thin- 
wire  probes/posts  inside  the  cavities  is  presented.  Coupled  integral  equations  are 
formulated  in  terms  of  the  electric  current  on  the  probes  and  the  equivalent  magnetic 
current  in  the  slots  and  a  numerical  solution  technique  based  on  the  moment  method  is 
used  to  solve  them. 

The  Ewald  method  is  employed  to  accelerate  the  convergence  of  the  free-space 
periodic  Green’s  function  appearing  in  the  kernel  of  the  integral  equations.  The  Ewald 
splitting  parameter  is  determined  using  a  special  method  (D.  Jackson,  private 
communication)  which  allows  for  good  accuracy  in  the  summations  over  a  wide  band  of 
frequencies.  The  need  to  further  speed  up  the  computation  led  to  the  usage  of  matrix 
interpolation  techniques  wherein  matrix  elements  are  computed  by  direct  means  at  only  a 
few  frequencies  and  then  interpolated  at  many  interior  frequencies.  Numerical  results  for 
the  input  impedance  of  the  probe  and  shielding  effectiveness  are  presented  for  various 
values  of  cavity  parameters. 

The  shielding  effectiveness,  a  means  of  determining  the  influence  of  the 
transmission  path  (and  environment)  on  the  shape  and  magnitude  of  a  transient  exciting 
signal,  is  evaluated.  It  is  hoped  that  this  study  will  shed  additional  light  on  how  the 
transmission  path  affects  the  characteristics  of  a  signal  which  enters  a  complex  electronic 
environment  and  reaches  an  interior  point  deep  inside  where  a  digital  circuit  may  be 


located. 
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Chapter  2 


AnalysisOf  Slot-Coupled  Cascaded 
Rectangular  Cavities 

Coupled  integral  equations  and  a  numerical  solution  technique  based  on  the 
moment  method  for  slot-coupled  cascaded  rectangular  cavities  are  presented  in  this 
chapter.  The  cascaded  rectangular  cavities  have  thin-wire  probes/posts  present  in  them. 
Free-space  periodic  Green's  functions  are  derived  for  the  equivalent  magnetic  current  in 
the  slot  and  the  electric  current  on  the  probe  for  a  single  rectangular  PEC  cavity.  The 
integral  equations  for  an  arbitrary  number  of  slot-coupled  cascaded  rectangular  cavities 
are  then  developed  as  a  general  case. 


1.  Introduction 

The  two  methods  that  are  commonly  used  for  electromagnetic  analysis  are  the 
mode-theory  formulation  [1,  2]  and  the  image-theory  formulation  [3],  with  the  latter 
being  employed  in  this  work.  The  free-space  periodic  Green's  function  appearing  in  the 
kernel  of  the  integral  equations  is  obtained  from  the  summation  of  the  elements  derived 
from  image-theory  analysis.  The  structure  of  primary  interest  is  a  series  of  slot-coupled 
cascaded  rectangular  cavities  with  an  arbitrary  number  of  sections  that  contain  thin-wire 
probes/posts  and  which  are  flanked  by  infinite  PEC  screens  at  both  ends.  In  general,  a 
plane-wave  or  a  current  excitation  may  be  used  to  excite  the  narrow  slots,  or  a  voltage 


7 


source  may  be  used  to  excite  the  probes  in  the  cavities.  The  analysis  of  this  structure 
requires  the  solution  of  coupled  integral  equations  which  include  the  free  space  periodic 
Green's  function  as  part  of  their  kernels.  The  coupled  integral  equations  are  solved  by  the 
method  of  moments  [4], 

2.  Analysis  Of  A  Single  Rectangular  Cavity 
To  analyze  the  structure  shown  in  Fig.  2-3  for  the  dipole  excitation,  one  must  first 
derive  the  Green's  function  for  an  jc  —  directed  magnetic  dipole  and  a  y-directed  electric 
dipole  in  a  rectangular  PEC  cavity  of  width  2a,  height  2 b ,  and  depth  c.  The  geometry 
for  this  structure  is  shown  in  Fig,  2-1. 


.v 


■fX 

■  2 a 

2b 

S l,  i.«d 

t'- 

Z 

c 

Fig  2-1.  An  x  —  directed  unity  strength  magnetic  dipole  and  a  y-directed  unity 
strength  electric  dipole  inside  a  rectangular  PEC  cavity. 

The  electric  vector  potential  Green's  function  for  the  x- directed  magnetic  dipole  located 
ata  point  (x  ',y',z')  inside  the  cavity  must  satisfy  the  differential  equation 

(V2  +  j32)Gx(x,y,z)  =  -e 8[x-x')8[y-y')8[z-z')  ,  (2,1) 


where 
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/3  =  UJy[/J£  (2.2) 

and  the  medium  inside  the  cavity  is  characterized  by  [ju.s] .  The  boundary  conditions  on 
the  cavity  walls  are 


— Gx(x,y,z)  =  0  for  x  =  ±a,  (2,3) 

By 


-^-Gx(x,y,z)  =  0  for  y  =  ±b  (2.4) 

oy 


and 


— Gx(x,  y,  z)  =  0  for  z  =  0and  z  =  c.  (2.5) 

oz 

In  the  image-theory  formulation  [3,  5]  the  source  is  first  imaged  across  the  cavity 
walls  as  shown  in  Fig.  2-2.  The  rectangle  shown  in  the  figure  with  the  dashed  lines 
represents  one  unit  cell  and  determines  the  periodicity  of  the  structure  in  the  *  and  y 
directions.  The  magnetic  dipole  shown  in  the  rectangle  with  solid  lines  is  the  original 
dipole  which  is  imaged  along  the  x  and  y  directions  repeatedly  to  obtain  the  periodicity 
seen  in  the  figure.  The  Green's  function,  which  satisfies  the  differential  equation  given  in 
(2.1),  the  boundary  conditions  given  in  (2.3),  (2.4)  and  (2.5),  and  the  radiation  condition, 
can  be  obtained  by  summing  the  free-space  Green's  function  for  each  dipole.  The  free 
space  periodic  Green's  function  for  an  x- directed  magnetic  dipole  located  at  a  point 
(x'.y'.z')  in  a  medium  characterized  by  [ju.s)  is  given  as 
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Thus,  if  all  the  contributions  are  counted,  the  image-theory  Green's  function  for  an 
x-directed  magnetic  dipole  inside  a  rectangular  PEC  cavity  is 
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where  the  8  contributions  are  a  result  of  8  magnetic  dipoles  in  one  unit  cell  as  shown  in 
Fig  2-4.  The  different  signs  on  the  contributions  are  due  to  the  x -directed  and  the  -x- 
directed  magnetic  dipoles,  as  shown  in  Fig.  2-2.  Also, 

R,nnpi  =J(£- mDx  f  +  (Vi  - nDy )2  +  ( c  -  pDz ) 2 ,  i  =  1,  ...  ,8  (2,9) 

and  the  terms  ,  7,. ,  and  £.  are 

4  =  &  =  £5  =4  =x-x',  (2.10) 

£2  =  £3  =4  =g7  =x+x'-2ci,  (2.11) 

m  =  rh=%=%  =  y-y' ■  (2'12) 

rh=ri4=rh  =%  =  y+ y'  -2b '  (2-13) 

Ci=C2=C3=^  =  z-z,l  (2.14) 

&=C6=Ci  =&=z  +  z'-2c,  (2.15) 

and  Dx,  Dy ,  and  D,  are  given  by 

Dx  =  4 a ,  Dv  =  8 b  and  D_  =  4c.  (2.16) 

The  terms  Dx,  Dv,and  Dz  represent  the  periodicity  of  the  structure  in  the  x,  y  and  z 
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directions,  respectively.  Similarly,  the  image-theory  Green's  function  for  a  y- directed 
electric  dipole  located  at  a  point  (x'.y'.z')  inside  a  rectangular  PEC  cavity  is  written  as 


G(r,r') 


E 

m,n,p=— oo 


g  2 

g  j  ^  ^rrinp  3 

1 

R  , 

mnp  1 

R  , 

mnp  2 

R  , 

mnp  6 

R 

mnp  4 

g  3  PR-mnp  5 

—jj3R  . 

C>  J '  mnp  o 

-j/3R  „  ' 

c>  J '  mnp  o 

1 

^ mnp  5 

R  „ 

mnp  D 

R  7 

mnp  7 

1 

R  * 

mnpo  j 

(2.17) 


where  the  terms  have  their  usual  meanings  and  the  8  contributions  are  a  result  of  8 
electric  dipoles  in  one  unit  cell  as  shown  in  Fig  2-4. 

In  Fig.  2-3(a),  the  structure  under  examination  is  a  rectangular  PEC  cavity 
between  two  infinite  PEC  screens  with  narrow  slots  in  them.  Inside  the  cavity  is  a  thin- 
wire  probe  which  may  connect  via  a  500  coaxial  cable  to  a  signal  source  having  an 
internal  impedance  of  500 .  The  side  view  of  the  structure  is  illustrated  in  Fig.  2-3(b).  As 
seen  from  the  figure,  the  length  of  the  thin-wire  probe  is  h  and  its  radius  is  r .  The  probe 
axis  is  located  in  the  yz  plane  with  d  representing  its  displacement  from  the  z  =  0 
plane.  The  slots  are  centered  on  each  face  of  the  rectangular  cavity  and  their  lengths  are 
2 Ly  and  2 L2,  respectively,  with  their  widths  being  2w.  The  width,  height,  and  depth  of 
the  cavity  shown  in  the  figure  is  2a,  2b  and  c,  respectively.  The  medium  inside,  in  the 
left  half-space  where  z  <  0 ,  and  in  the  right  half-space  where  z  >  c ,  is  characterized  by 
(//,£•)  . 

The  complete  3D  unit  cell  due  to  imaging  is  shown  in  Fig,  2-4.  In  the  figure,  the 
cavity  shown  with  solid  lines  is  the  original  rectangular  cavity  and  the  cavities  shown 
with  dashed  lines  represent  the  image  cavities.  The  electric  current  on  the  probe  and  the 
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magnetic  current  in  the  slot  are  imaged  simultaneously.  The  resulting  unit  cell  has  8 
magnetic  dipoles,  as  shown  in  Fig.  2-4,  due  to  the  images  of  slot  current  Mvl  on  the 

cavity  wall.  The  images  of  the  slot  current  on  cavity  wall  2,  designated  as  mx2 , 
contribute  to  4  magnetic  dipoles  in  the  unit  cell.  The  current  ly  on  the  probe  has  8 

contributions  in  the  unit  cell.  A  rectangular  cavity  with  2  narrow  slots  is  pictured  in  Fig. 
2-3  to  show  the  presence  of  8  magnetic  dipoles  due  to  wall  1  and  4  due  to  wall  2  in  the 
same  unit  cell.  The  Green's  function  for  an  x-  directed  magnetic  dipole  on  slot  2  inside  a 
rectangular  PEC  cavity  is  presented  in  the  next  section. 

The  xy  and  xz  plane  imaging  is  shown  in  Fig.  2-5  for  a  single  cavity  which  can 
be  generalized  to  an  arbitrary  number  of  cascaded  rectangular  cavities.  The  rectangle  in 
darker  solid  lines  represents  the  original  cavity  and  it  is  imaged  about  the  xy  and  the  xz 


planes  successively. 
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(b) 

Fig.  2-3.  A  rectangular  PEC  cavity  between  two  infinite  PEC  screens  with  narrow 

slots  in  them  and  with  a  thin-wire  probe  inside,  (a)  perspective  view  (b)  side 
view 
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Fig.  2-5.  Equivalent  model  for  interior  region  with  imaging  about  xy  and  xz  planes 
(single  rectangular  PEC  cavity  with  two  slots  and  one  probe)  -  side  view 
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3.  Arbitrary  Number  Of  Slot-Coupled  Cascaded  Rectangular  Cavities 

A  general  case  of  an  arbitrary  number  /  of  slot-coupled  cascaded  rectangular 
cavities  is  pictured  in  Fig.  2-6.  Each  cavity  has  a  thin-wire  probe/post  in  it  and  any  probe 
can  be  connected  via  a  50Q  coaxial  cable  to  a  signal  source  having  an  internal 
impedance  of  500 .  The  last  rectangular  section  is  the  J'h  cavity  and  the  sections 
between  the  first  and  the  last  rectangular  cavity  are  denoted  by  j .  The  width,  height,  and 
depth  of  each  cavity  shown  in  the  figure  is  2a,  2b  and  c,  respectively.  The  length  of  the 
/h  thin-wire  probe  is  hj  and  its  radius  is  r .  The  probe  axes  are  located  in  the  yz  plane 

with  d1,  ... ,  (j -1)* c  +  dj ,  ...,  U-l)*c+  dj  representing  their  displacements  from  the 

z  =  0  plane.  The  length  of  the /h  slot  is  2Ly  and  its  width  is  2 w.  The  medium  inside  each 
cavity,  in  the  left  half-space  where  z< 0,  and  in  the  right  half-space  where  z>(J*c)  is 
characterized  by  (//,£•) . 

The  procedure  followed  in  the  previous  section  is  utilized  again  and  the  general 
problem  of  an  arbitrary  number  of  slot-coupled  cascaded  rectangular  cavities  is  broken 
down  into  single  cavities.  The  treatment  is  essentially  the  same  and  the  unit  cell  obtained 
due  to  image  theory  in  Fig.  2-4  is  used  separately  for  each  cavity.  The  procedure  by 
which  one  may  derive  the  coupled  integral  equations  for  a  thin-wire  behind  an  infinite 
PEC  screen  is  discussed  in  [5].  In  [3],  Kustepeli  describes  the  coupled  integral  equations 
for  a  narrow  slot  in  an  infinite  PEC  screen  backed  by  a  single  rectangular  PEC  cavity 
with  a  thin-wire  probe  inside  it. 
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The  formulation  discussed  in  [3]  has  been  extended  in  this  work  to  an  arbitrary 
number  of  slot-coupled  cascaded  rectangular  cavities.  As  an  aid  in  the  formulation  of  the 
coupled  integral  equations  for  the  structure  of  interest,  the  equivalent  models  for  the 
interior  and  exterior  regions  of  the  structure  are  developed.  The  tangential  magnetic  field 
is  then  forced  to  be  continuous  in  the  apertures  and  the  total  tangential  electric  field  on 
thethin-wire  probes  in  each  cavity  is  made  zero.  The  resulting  coupled  integral  equations 
derived  for  the  electric  current  on  the  probes  and  the  equivalent  magnetic  current  in  the 
slots  are  presented  for  the  structure  shown  in  Fig.  2-6,  The  equations  (2.18)  and  (2.21)  are 
the  coupled  integral  equations  for  slot  1  and  slot  7+1  where  the  slot  has  been  shorted 
and  a  short-circuit  magnetic  current  is  placed  in  its  position.  Equation  (2.20)  is  the 
coupled  integral  equation  for  slot  j  where  1<  j<  7+1.  The  total  tangential  electric  field 
on  thethin-wire  probes  is  made  zero  and  the  resulting  coupled  integral  equation  is  listed 
in  (2.19). 


I — 1 

[  d 2 

jlirujii 

| 

dx2 

\ 

+  p2  J  ml(x,){K'i(x-x,)  +  Ksc^(x,x,)}dx' 


x'=—Li 


+  Y~  f  I1{y')G?{x,y')dy' 


47T  / 

y  =-h 


i — 1 

d2 

j  2tt uif-i 

| 

dx2 

V 

J  m2(x')KsJti[x,x')dx' 


x'——L) 


=  -H[c ,  on  slot  1. 
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(2,20) 


(2.21) 
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In  equations  (2.18-2,21),  the  axial  variation  of  the  magnetic  current  on  slot  j  is  given  by 
nij  ( x)  and  the  electric  current  on  the  thin-wire  probe  j  is  given  by  lj(y) .  In  equation 

(2.19) ,  the  voltage  at  the  center  of  probe  j  is  given  by  Vr  The  delta  function,  6(y+b)  in 

(2.19)  specifies  that  Vl  on  probe  j  is  located  at  the  y  =  -b  plane.  The  short-circuit 
magnetic  field,  in  (2.18),  (2.20)  and  (2,21)  is  the  field  at  the  location  of  slot  j  when 
it  is  shorted. 

The  contribution  to  the  kernel  due  to  slot  1,  K[  in  (2.18)  is  written  as 

Ai 


<  (jc-jcO 


(2,22) 


where 


«/  =  J(W)!  +  4(^)!sin2y  (2.23) 

In  equation  (2.21),  K'Sj  is  the  kernel  due  to  slot/+l  and  it  can  be  defined  like  K's  in 
(2,22).  The  third  term  of  (2.20),  TQ  due  to  the  presence  of  slot  j  in  cavity  j  is  written 

cvj 

as 


KSj  =  [ 

cvi  J 


II 
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J  I  2  a 
y'=-wyW  -y 


{GSi  (x,0,z',x',  y',z') 


+GSl(x,0,  z',x',-2b  +  y',z'))dy',  j  =  1,  .  ..,7+1. 


(2.24) 


The  Green's  function  GSJ  in  (2.19)  is  given  by 
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G'd  ( y,  x\  y')  =  ^-Gs‘  (0,  y,  z\  x' ,  y' ,  z') 
oz 

+  -^-G'J(0,  y,z)x',-2b  +  y’.z'),  7=1,  ....  J, 

oz 


(2.25) 


where  GSi  in  (2.24)  and  (2,25)  is  equivalent  to  Gx  in  (2.8).  In  (2,20)  K’J?  due  to  the 
presence  of  slot j+l  in  cavity  j  is  written  as 
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and  G'li+1  in  (2.19)  is  written  as 
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In  (2,26)  and  (2,27),  GSi+1  is  the  image-theory  Green's  function  for  an  x  —  directed 
magnetic  dipole  inside  a  rectangular  PEC  cavity  and  is  written  as 


GSi+t 

where  the  4  contributions  are  as  a  result  of  4  magnetic  dipoles  in  one  unit  cell  as  shown  in 
Fig  2-4  for  the  magnetic  current  on  slot  2  and 
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Knpib  =  J  (L  -  mDx  )2  +  (Vib  -  nDy  f  +  (Cb  -PD:)2  ’ 


i  =1,2, 3, 4.  (2,29) 


21 


where  D.  D,  ,  and  D,  are  defined  in  (2.16)  and 

X  (/  z 


4  =4,  =  *-*'. 

(2.30) 

4ft  =  4,  =  *  +  *'- 2a, 

(2.31) 

dv,  =  dib  =  y-  y'' 

(2.32) 

rhb  =  rub  =  y+y'-2b' 

(2.33) 

C 2b  ^3b  ^ ^ b  ^  " 

(2.34) 

The  kernel  due  to  the  presence  of  the  probe  inside  the  cavity,  Kpj  ,  is  written  as 
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and  Gp  in  (2,18)  and  (2.20)  is  written  as 
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In  (2.35)  and  (2.36),  Gp>  is  the  image-theory  Green's  function  for  a  y- directed  electric 
dipole  inside  a  rectangular  PEC  cavity  and  is  similar  to  Gy  in  (2.17). 


4.  Numerical  Solution 

The  coupled  integral  equations  given  in  (2.18-2.21)  are  solved  by  the  method  of 
moments  [4],  In  this  method  the  unknown  magnetic  currents  m} (x)  and  the  unknown 

electric  currents  /.(y)  are  expanded  in  terms  of  the  pulse  basis  functions,  each  weighted 

by  unknown  coefficients  m  and  i  ,  respectively.  The  triangle  functions  are  used  for 

testing.  Detailed  information  about  expansion  and  testing  procedures  is  given  in  [4,  7], 
The  resulting  matrix  equation  is  written  as 
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where  Ymnj,  Ymnu+i,  and  Ymnj+u  are  the  admittance  terms  due  to  the  slots.  The  elements  of 

the  sub-matrices  representing  coupling  from  the  probe  to  the  narrow  slot  are  given  by 
C™  and  C™1  represents  the  slot-to-probe  coupling  terms.  The  impedance  sub-matrix 

term  for  the  probe  is  Z'2 .  The  expression  for  the  elements  of  the  admittance  sub-matrix 
Ymni  is  given  by 


Y  =  Ys> 

mri:  mri: 


+  r\  i=l  or  7, 

nin  ’  • 


(2.38) 


where 
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In  (2.39)  xm  and  x„  can  be  determined  from 


xp=-Li  +  pASj,  p  =  l,2,...,N 


(2.40) 


with  p  standing  for  m  or  n .  The  sub-domain  segment  length  As  is  given  by 


2A 


(N  +  l) 


(2.41) 


In  (2.40)  and  (2.41),  N  is  the  number  of  basis  functions  used  to  approximate  the 
magnetic  current  in  the  slot.  One  can  write  Yr2  as 
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The  expression  for  the  elements  of  the  admittance  sub-matrix  Fm„  is  given  by 


F  =  fcv'  z=l  / 

mn/+l,z  mn/+l  '  '  ■■■'  U  ’ 


(2.44) 


where  F";  and  elements  of  the  sub-matrix  Ymn  appearing  in  (2.37)  are  defined  like 


C  of  (2.42). 


In  (2.37),  the  elements  of  the  sub-matrix  C^f  representing  the  coupling  from  the 
probe  to  the  slot  is  given  by 
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The  width  of  a  pulse  basis  function  used  to  expand  the  electric  current  on  the  probe  is 


def i  ned  as 
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2 h. 


(P  +  1) 


(2.47) 


where  P  in  (2.46)  and  (2.47)  is  the  number  of  unknowns  used  to  expand  the  current.  The 
elements  of  the  sub-matrix  Csmf  '‘  are  defined  like  C^f  in  (2,45). 

In  (2.37),  the  elements  of  the  impedance  sub-matrix  Zpmn  is  written  as 
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where 
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The  expression  for  the  elements  of  the  sub-matrix  C^‘  appearing  in  (2.37) 
representing  the  coupling  from  the  slot  to  the  probe  is  given  by 
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(2,50) 


The  elements  of  the  sub- matrix  C,^i+1  are  defined  like  C™  in  (2.50). 

Finally,  the  elements  of  the  excitation  vector  in  (2.37)  are  given  by 


VL  = 


-V,  m  =  P*  [i  +  J)  + 


P  +  1 


(2.51) 


0, 


otherwise. 
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For  plane-wave  excitation  of  the  slot,  the  term  lm  of  (2.37)  is  given  by 

Imi  =  A sHs'[xm),  m  =  (i-l)*P  +  (1 . P)  (2.52) 

The  approximation  [5]  employed  in  (2.52)  is  a  very  good  one,  provided  Hsxc  does  not 
vary  rapidly  over  an  interval  of  2As  . 


5.  Conclusions 

In  this  chapter  coupled  integral  equations  and  a  numerical  solution  technique 
based  on  the  moment  method  for  slot-coupled  cascaded  rectangular  cavities  are 
presented.  Free-space  periodic  Green's  functions  for  a  magnetic  dipole  and  an  electric 
dipole  in  a  rectangular  PEC  cavity  are  derived  and  coupled  integral  equations  for  an 
arbitrary  number  of  slot-coupled  cascaded  rectangular  cavities  are  then  presented  as  a 
general  case.  In  this  chapter,  a  general  solution  is  described  as  a  precursor  to  a  numerical 
solution  for  single,  double,  and  triply  cascaded  rectangular  cavities  which  are  discussed 
in  detail  in  Chapter  5. 
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Chapter  3 


Ewald  Method  And  Determination  Of  Splitting 

Parameter 

The  free-space  periodic  Green's  function  appears  as  the  kernel  of  the  integral 
equations  formulated  for  slot-coupled  cascaded  rectangular  cavities.  The  Ewald  method 
[1,  2]  can  be  used  to  efficiently  accelerate  the  convergence  of  the  infinite  series 
encountered  in  the  calculation  of  matrix  elements  obtained  by  employing  the  method  of 
moments  for  the  numerical  solution  of  the  integral  equations.  In  this  method,  the  free- 
space  periodic  Green's  function  is  expressed  as  the  sum  of  "spectral"  and  "spatial"  series. 
Each  of  these  series  possesses  Gaussian  decay,  leading  to  a  series  representation  that 
typically  has  a  very  fast  convergence  rate.  The  spectral  and  spatial  series  are  written  in 
terms  of  the  complementary  error  function,  involving  a  splitting  parameter E.  The 
method  of  determining  the  splitting  parameter  and  numerical  results  for  the  different 
choices  of  the  splitting  parameter  are  presented  in  this  chapter. 

1.  Introduction 

The  Ewald  method  was  first  developed  and  used  in  1921  by  P.P.  Ewald  [1]  to 
evaluate  the  summations  occurring  in  crystal  lattice  theory.  In  the  evaluation  of  the  free- 
space  periodic  Green's  function  for  two-dimensional  and  three-dimensional  cases  with 
the  Ewald  method,  the  spectral  and  the  spatial  series,  which  are  written  in  terms  of  the 
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complementary  error  function,  involve  a  splitting  parameter E.  In  [9]  it  is  stated  that  E  is 
an  arbitrary  number  with  an  optimum  value,  Eopt,  given  to  balance  the  convergence  rate 

between  these  two  series,  thereby  minimizing  the  total  number  of  terms  needed  for  the 
calculation.  If  E  is  increased  beyond  Eopt ,  then  successive  terms  in  the  spatial  series 

decay  faster  while  successive  terms  in  the  spectral  series  decay  more  slowly.  In  studies 
involving  the  Ewald  method  [9,  11],  the  Eopt  value  is  generally  used  to  sum  the  spectral 

and  the  spatial  series  arising  in  the  analysis  of  structures  having  small  periodic  spacing. 
Although  it  is  stated  in  [9]  that  the  splitting  parameter  is  arbitrary,  in  [8],  Kustepeli  and 
M  artin  have  shown  that  one  cannot  get  correct  results  even  when  E  =  Eopt  for  structures 

having  large  periodic  spacing.  Therefore,  the  splitting  parameter  should  always  be  chosen 
carefully. 

One  can  apply  the  Ewald  method  in  calculating  the  free-space  periodic  Green's 
function  due  to  three-dimensional  periodic  arrays  [1].  For  the  three-dimensional  case  in 
[10],  the  value  of  E  used  in  the  computations  is  determined  by  examining  only  the 
convergence  of  the  spatial  series.  For  two-dimensional  arrays,  the  choice  of  this 
parameter  is  very  important  as  in  the  three-dimensional  case,  to  balance  the  convergence 
rate  in  the  analysis  of  structures  with  small  periodic  spacing  and  to  obtain  correct  values 
in  the  analysis  of  structures  having  large  periodic  spacing.  One  can  also  use  the  Ewald 
method  for  the  three-dimensional  case  after  applying  image- theory  to  the  structure. 

2.  Determination  of  Splitting  Parameter 
In  [8],  Kustepeli  and  M  artin  present  an  expression  for  Eopl  which  can  be  used  in 

the  calculation  of  free-space  periodic  Green's  function  for  the  three-dimensional  case 
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having  small  periodic  spacing.  Results  obtained  by  the  Ewald  method  are  compared  with 
those  obtained  by  the  well-known  Shanks'  transform  [4,  5]  computed  using  Wynn's  s- 
algorithm  [6].  One  of  the  most  effective  means  for  evaluating  the  Green's  function  given 
in  (2.6)  of  the  previous  chapter  is  the  Ewald  Method  [1],  In  this  method,  the  Green's 
function  is  expressed  as  the  summation  of  two  series  such  that 

G  =  Gj  +  G2  .  (3.1) 

In  (3.1),  Gj  and  G2  utilize  the  complementary  error  function,  erfc[z) .  By  virtue  of  the 
presence  of  the  complementary  error  function  in  these  series,  a  very  rapid  convergence 
rate,  resulting  from  the  Gaussian  decay  found  in  the  terms,  is  achieved.  Since  Gaussian 
decay  is  involved  in  each,  the  value  of  the  summation  is  generally  obtained  in  a  very 
small  number  of  terms  if  E  is  properly  chosen.  For  the  three-dimensional  case,  the 
Green's  function  in  (2.6)  is  expressed  as  a  sum  of  two  series  as  given  in  (3.1).  The 
spectral  series  G1  is  written  as  [3,  8] 
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and  Dx ,  Dy  and  Dz  are  given  in  (2.16).  The  spatial  series  G2  is  given  by  [3,  8]  as 
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which  can  also  be  written  as 
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where  is  defined  in  (2.9).  The  optimum  splitting  parameter  [3,  8]  used  to  balance 
the  convergence  rates  of  Gl  and  G2  for  the  three-dimensional  (3D)  array  is 
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The  complementary  error  function  appearing  in  the  equations  (3.4)  and  (3.5)  is 
given  by  [3,  8] 


erfc(z)  =1  -erf(z)  = 


e  '  dt 


where  erf(-),  the  error  function,  is  given  by 


erf(z)  = 


(3.7) 


(3.8) 


By  using  the  following  relation, 
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erfc(z)  =  e  z 2  w(jz)  ,  (3.9) 

the  complementary  error  function  appearing  in  the  equations  (3.4)  and  (3.5)  is  calculated. 
The  function  w(  •)  in  (3.9)  is  utilized  to  compute  the  complementary  error  function  for  a 
complex  argument.  This  function,  which  was  first  devised  and  tabulated  by  Faddeyeva 
and  Terent'ev  [12],  can  be  computed  using  a  code  written  by  Poppe  and  Wijers  [13, 14], 
The  method  explained  in  [8]  to  determine  El Eopt  involves  a  splitting  parameter 

which  is  said  to  be  dependent  on  periodic  spacing  of  a  structure  alone  and  the  choice  of 
which  is  more  or  less  arbitrary.  Jackson  suggested  a  method  [7]  wherein  the  splitting 
parameter  is  not  treated  as  arbitrary  and  is  said  to  be  dependent  on  periodic  spacing,  the 
number  of  significant  figures  lost  in  the  calculation  of  the  Green's  function  and  the 
frequency. 

The  ''optimum"  E  parameter,  Eopt ,  that  results  in  the  same  asymptotic  rate  of 
decay  for  Gl  and  G2  is  given  in  (3.6).  At  high  frequencies,  one  gets  very  large  nearly- 
equal-in-magnitude  numbers  of  opposite  signs  for  the  first  several  terms  of  Gl  and  G2  . 
These  terms  largely  cancel  one  another,  resulting  in  a  sum  that  is  of  moderate  value. 
However,  this  cancellation  generally  results  in  an  apparent  convergence  to  incorrect 
values  and  a  loss  of  significant  figures.  To  avoid  this  problem,  one  needs  to  limit  the  size 
of  the  largest  terms  in  both  series.  Jackson  [7]  devised  a  method  for  choosing  E  that  is 
larger  than  the  optimum,  which  solves  this  problem. 

The  goal  is  to  limit  the  size  of  the  largest  term  relative  to  that  of  the  overall 
Green's  function.  The  largest  terms  in  each  series  come  from  the  (0,0,0)  terms.  One 


chooses  a  limit  of 
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(3.10) 


and 


|G2i00o  |  <  10"  |G|,  (3.11) 

where  G:  and  G2  have  their  usual  meanings,  the  subscript  (0,0,0)  refers  to  the  first  term 
of  the  series  (m  =  n  =  p  =  0)  and  the  parameter " d "  indicates  the  significant  figures  lost 
in  the  calculation  of  the  Green's  function.  The  magnitude  of  the  overall  Green's  function 
is  approximated  as 
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This  approximation  is  a  reasonable  one  unless  the  distance  between  the  source  and 
observation  points  becomes  comparable  to  the  distance  between  the  source  point  and  the 
boundary  of  the  unit  cell  (so  that  image  terms  are  important).  Even  then,  the  above 
estimate  is  probably  correct  to  within  an  order  of  magnitude. 

In  Jackson's  method  one  uses  the  asymptotic  form  of  the  complementary  error 
function  (valid  at  high  frequencies)  and  using  (3.10)  writes  Gl  of  (3.2)  as 
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One  can  define 
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then  one  gets 
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An  accurate  solution  to  the  transcendental  equation  in  (3.15)  is  given  by 
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The  inequality  in  (3.13)  has  to  be  maintained  and  hence  the  restriction  that 


E>-°^,  (3.18) 

-T. 

rendering  (3.13)  true. 

The  analysis  for  the  G2  series  is  similar  using  Jackson's  method.  One  uses  the 
asymptotic  form  of  the  complementary  error  function  (valid  at  high  frequencies)  and 


writes 
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One  can  write 
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then  (3.19)  becomes 
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The  solution  to  the  transcendental  equation  in  (3.21)  is  given  by 


x2=  ln(c2)  +  ln  J(c2)  +  ln  J(c2)  +  \n[Mc. 


(3.23) 


The  inequality  in  (3.19)  has  to  be  maintained  and  hence  the  restriction  that 


(3.24) 
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Jackson's  method  thus  chooses  an  overall  E  by  comparing  (3.6),  (3.18)  and  (3.24)  and 
writing 


f 
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K  ustepeli  and  M  artin  [8]  and  J  ackson  [7]  have  shown  that  E  =  Eopt  may  not  be  sufficient 

for  large  periodic  spacing  and  E  has  to  be  increased  beyond  this  value  to  achieve 
accurate  results.  The  reason  E  —  Eopt  does  not  work  for  larger  periodic  spacing  can  be 

explained  as  follows.  For  large  arguments  the  complementary  error  function  behaves  as 
exp(-z2)/z  which  is  the  reason  for  using  this  method  for  the  calculation  of  the  free- 
space  periodic  Green's  function.  For  large  periodic  spacing,  Eopt  given  by  (3.6)  becomes 
small,  and  the  imaginary  part  of  the  argument  of  erfc(-)  becomes  large  and  dominantfor 
the  first  several  terms  of  the  series.  Asa  consequence,  erfc(  •)  becomes  very  large  and, 
therefore,  one  gets  very  large  numbers  for  the  first  several  terms  of  these  series.  Similar 
comments  can  be  made  for  the  series  G:  of  (3.2)  that  includes  the  Gaussian  function 
directly.  Gaussian  decay  is  again  achieved  for  terms  with  large  indices,  and  the  two  series 
G1  and  G2  converge  to  very  large  nearly-equal-in-magnitude  numbers  of  opposite  signs. 
Since  the  values  of  the  spectral  and  the  spatial  sums  are  very,  very  large  but  have 
opposite  signs,  accuracy  is  seriously  compromised  when  adding  the  sums  of  the  two 
series  because  of  finite  precision.  The  result  is  the  apparent  convergence  to  incorrect 
values.  By  increasing  E  beyond  Eopt  by  a  sufficient  amount,  one  obtains  small  values 

for  the  imaginary  part  of  the  argument  of  the  complementary  error  function  and  small 
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positive  values  for  -a2mnpl  E2  in  (3.2)  for  the  first  several  terms  of  the  series.  Asa 

result,  one  avoids  the  inaccuracy  resulting  from  adding  numbers  of  nearly  equal 
magnitude  but  of  opposite  sign,  and  a  correct  sum  is  obtained  for  G .  Thus  a  value  for  E 
is  chosen  that  is  the  maximum  of  (3.6),  (3.18)  and  (3.24). 

3.  Results 

The  following  section  compares  Jackson's  method  [7]  of  determining  E  with  an 
arbitrary  method  of  determining  E  and  presents  comparative  results  for  the  same.  The 
geometry  chosen  for  this  analysis  is  shown  in  Fig.  3-1.  As  shown  in  the  Fig.  3-l(a),  the 
structure  under  examination  is  a  narrow  slot  in  an  infinite  PEC  screen  backed  by  a 
rectangular  cavity  with  a  thin-wire  probe  present  inside  the  cavity.  The  side  view  of  this 
structure  is  illustrated  in  Fig,  3-l(b).  The  length  of  the  thin-wire  probe  is  h  with  its  radius 
specified  as  r.  The  probe  axis  is  located  in  the  yz  plane  with  d  representing  its 
displacement  from  the  z  =  0  plane.  The  slot  length  is  2 L  and  its  width  is  2 w.  The 
width,  height,  and  depth  of  the  cavity  shown  in  the  figure  are  2a,  2b  and  c ,  respectively. 
The  medium  inside  the  cavity  and  in  the  left  half-space  where  z<0  is  characterized  by 
(ju,e).  The  probe  is  connected  via  a  50fi  coaxial  cable  to  a  signal  source  having  an 
internal  impedance  of  50fi .  The  structure  parameters  are  shown  in  T able  3.1. 

The  real  and  imaginary  parts  of  the  input  impedance  of  the  thin-wire  probe  for  the 
structure  given  in  Table  3.1  in  the  frequency  range  between  6  GHz  and  12  GHz  are 
shown  in  Fig.  3-2,  In  this  figure,  the  measured  data  from  [3]  is  compared  with  the 
numerical  results  obtained  by  using  Jackson's  method  and  E  =  E  .  The  figure  shows 

that  the  numerical  results  begin  to  fail  when  the  frequency  is  increased  while  E  =  Eopl . 
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As  shown  in  the  figure,  the  real  part  of  the  computed  input  impedance  becomes  even 
negative  for  some  frequencies  when  E  =  E  ;  in  addition,  the  results  for  the  imaginary 

part  are  not  even  in  the  range  of  the  figure  above  10.4  GHz,  showing  a  large  error  in  the 
computation  of  these  values.  Since  the  value  of  the  splitting  parameter  is  not  appropriate 
for  the  calculation  of  the  matrix  elements,  error  occurs  in  the  calculation  of  the  probe 
current  and  input  impedance.  The  correct  values  of  the  matrix  elements  can  be  obtained 
by  increasing  the  splitting  parameter  by  a  sufficient  amount,  which  results  in  the  correct 
values  for  the  probe  current  and  hence  the  input  impedance.  As  shown  in  the  figure,  by 
using  the  J  ackson's  method,  one  can  get  very  good  agreement  between  the  measured  [3] 
and  the  numerical  results.  Thus  one  can  conclude  from  Fig  3-2,  the  optimum  splitting 
parameter  ( E  =  E )  may  not  be  suitable  for  some  structures.  Therefore,  the  splitting 

parameter  should  be  chosen  carefully  and  one  may  use  Jackson's  method  to  determine 
the  E  parameter  without  the  need  for  guess  work  or  trail  and  error. 

Variation  of  E/Eopl  with  J  ackson's  E  is  smooth  in  the  entire  frequency  range 

as  shown  in  Fig.  3-3  compared  to  the  arbitrary  usage  of  E  =  E  t  and  E  =  2*  E  .The 

other  major  consideration  for  employing  Jackson's  method  was  CPU  time  reduction. 
CPU  times  for  numerical  results  obtained  by  using  E  —  E,  and  E  =  2*E  ,  were 

compared  to  results  obtained  by  using  Jackson's  E .  As  shown  in  Fig.  3-4,  cumulative 
computation  time  is  least  for  E  =  E  and  it  increases  as  the  multiplication  factor 

increases.  Also,  cumulative  CPU  time  for  Jackson's  method  is  between  E  =  E  and 

E  =  2*  Eopt  in  the  frequency  range  of  interest.  Thus,  J  ackson's  method  is  best  suited  in 

terms  of  CPU  time  and  accuracy  of  results. 
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Table  3.1 

Structure  Parameters 


Parameter 

Structure  A 

Cavity  width,  2 a 

22.86 

Cavity  height,  2 b 

10.16 

Cavity  depth,  c 

41.275 

Slot  length,  2 L 

15.0 

Slot  width,  2 w 

0.5 

Probe  length,  h 

6.35 

Probe  radius,  r 

0.254 

Probe  position,  d 

17.4625 

All  parameters  are  in  millimeters. 


TABLE3.2 

Resonant  Frequencies" 


Structure  A 

TEm 

7.50 

TE102 

9.79 

TEW3 

12.72 

TE2Q1 

13.61 

TEm 

15.94 

Probe  half-wave  resonance 

11.81 

Slot  half-wave  resonance 

10.00 

All  frequencies  are  in  GHz. 


Structure  A 


A  narrow  slot  in  an  infinite  PEC  screen  backed  by  a  rectangular  cavity 
with  athin-wire  probe  in  it  -  (a)  perspective  view,  (b)  side  view. 
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(a) 


(b) 


Fig.  3-2  -  Input  impedance  at  the  base  of  the  probe  for  Structure  A  (a)  Real  part  of 
impedance,  (b)  Imaginary  part  of  impedance. 
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Frequency  (GHz) 


. E=Eopt 

*  Jackson's  E 
- E=2*Eopt 


Fig.  3-3.  Variation  of  E/Eopt  with  frequency  for  structure  A. 


Fig.  3-4. 


Cumulative  CPU  time  vs.  frequency  in  computing  input  impedance  at  the 
base  of  the  probe  for  structure  A . 
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4.  Conclusions 

In  this  chapter  a  discussion  of  the  Ewald  method  [1,  2]  which  is  used  to  efficiently 
accelerate  the  convergence  of  the  infinite  series  encountered  in  the  calculation  of  matrix 
elements  obtained  by  employing  the  method  of  moments  for  the  numerical  solution  of  the 
integral  equations  is  discussed.  Also,  Jackson's  method  of  determining  the  splitting 
parameter  in  the  Ewald  method  has  been  discussed  and  results  showing  the  accuracy  and 
efficiency  of  this  approach  compared  to  an  arbitrary  method  of  choosing  E  have  been 
detailed.  Expressions  for  choosing  the  correct  E  parameter  are  arrived  at  after  going 
through  formulations  from  Jackson's  method.  The  results  obtained  by  using  Jackson's 
method  are  compared  to  the  measured  results  [3]  and  accuracy  is  established. 
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Chapter  4 


Efficient  Modeeing  Using  Matrix 
Interpolation  Methods 

The  free-space  periodic  Green's  function  appears  as  the  kernel  of  the  integral 
equations  formulated  for  slot-coupled  cascaded  rectangular  cavities.  One  usually 
encounters  extremely  slowly  convergent  series  in  the  calculation  of  the  matrix  elements 
occurring  in  the  method  of  moments  (MoM)  numerical  solution  of  these  integral 
equations.  One  is  thus  left  with  a  computationally  heavy  task  of  computing  the  matrix 
elements.  One  approach  that  is  useful  for  this  task  is  to  employ  matrix  interpolation 
techniques  (M  ITs)  to  compute  the  elements  of  the  system  matrix  to  reduce  the  overall 
computation  time.  In  this  chapter  a  discussion  of  the  M  IT  and,  in  particular,  the  ratio  of 
polynomials  (ROPs)  method  [1]  is  provided  and  numerical  results  are  given  to  establish 
the  accuracy  and  efficiency  of  this  approach  compared  to  direct  computations. 

1.  Introduction 

An  investigation  of  an  efficient  method  for  the  computation  of  the  broadband 
performance  of  slot-coupled  cascaded  rectangular  cavities  using  frequency  interpolation 
[1,  3]  of  the  MoM  matrix  is  presented.  In  the  evaluation  of  cascaded  cavities  for 
obtaining  wideband  data,  the  MoM  numerical  solution  can  take  a  very,  very  long  time 
since  the  system  matrix  has  to  be  computed  separately  for  each  new  frequency.  With  an 
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interpolation  method,  the  elements  of  only  a  few  matrices  at  relatively  large  frequency 
intervals  are  directly  computed.  These  matrices  are  then  used  to  interpolate  the  elements 
of  the  system  matrix  at  many  intermediate  frequencies.  Interpolation  reduces  the  time  [1, 
3]  required  to  compute  the  cavity  performance  over  a  wide  frequency  band, 

T wo  methods  that  are  commonly  used  for  electromagnetic  analysis  are  the  finite- 
difference  time-domain  (FDTD)  method  [4,  5]  and  the  method  of  moments  (MoM  )  [6,  7, 
8],  FDTD  predicts  performance  over  a  wide  band  of  frequencies  in  one  simulation.  The 
frequency-domain  M  oM  approach  computes  the  system  matrix  one  frequency  at  a  time 
and  requires  the  computation  of  many  frequency  points  for  broadband  performance 
evaluation.  This  can  take  a  long  time  in  situations  where  each  matrix  element  requires  a 
large  amount  of  CPU  time  since  the  system  matrix  will  have  to  be  recomputed  for  each 
new  frequency  point.  Popular  methods  of  reducing  the  computation  time  and  significantly 
speeding  up  the  overall  simulation  process  rely  on  the  use  of  numerical  and  geometrical 
approximations  to  quickly  fill  the  system  matrix  or  efficient  matrix  inversion  and  solution 
algorithms  [7,  8], 

The  [Z]  interpolation  method  was  originally  proposed  by  Newman  and  Forrai  [9] 
for  the  scattering  analysis  of  a  microstrip  patch.  Newman  [3]  also  used  it  for  the 
impedance  analysis  of  a  straight  dipole  antenna  and  flat  square  plate.  In  [1,  2]  Virga  and 
Rahmat-Samii  discuss  [Z]  and  [Y]  interpolation  methods  that  involve  simple 
interpolation  functions  which  require  only  a  few  coefficients  and  which  can  be  applied  to 
a  wide  variety  of  geometrical  configurations.  The  technique  is  also  independent  of 
excitation  and  imposes  few  restrictions  on  geometry.  Their  work  expanded  the  utilization 
of  [Z]  interpolation  to  the  wideband  performance  evaluation  of  complex  antennas  and 
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also  investigates  a  comparable  [F]  interpolation  approach.  The  methods  are  applied  in 
the  context  of  a  triangular  surface  patch  M  oM  formulation  [2],  In  [1],  the  [Z]  and  [Y] 
interpolation  techniques  are  compared  and  the  accuracy  and  computational  efficiency  of 
the  two  methods  are  discussed  along  with  the  implementation  guidelines. 

Generally,  the  elements  of  [Z]  and  [Y]  calculated  in  free-space  vary  slowly  with 
frequency  while  the  elements  calculated  in  the  presence  of  the  rectangular  cavity 
fluctuate  rapidly  with  frequency.  The  elements  of  [Z]  and  [F]  calculated  in  the  presence 
of  the  rectangular  cavity  strongly  depend  upon  the  overall  behavior  of  the  entire  structure. 
The  elements  of  [Z]  and  [F]  in  free-space  can  be  evaluated  over  a  frequency  range  by 
simple  and  low-order  interpolation  functions,  such  as  a  quadratic  polynomial.  In  contrast, 
the  interpolation  of  elements  in  the  presence  of  the  rectangular  cavity  over  a  frequency 
range  requires  interpolation  functions  such  as  the  ratio  of  two  nth  order  polynomials  to 
accurately  capture  the  frequency  behavior.  The  interpolation  method  used  to  compute  the 
system  matrix  discussed  in  Chapter  2  of  this  report  is  outlined  in  the  next  section. 

2.  Methodology  -  Ratio  of  Polynomials  Interpolation 
The  ROPs  interpolation  method  utilized  here  focuses  on  matrix  interpolation 
methods  that: 

.  are  easily  implemented  in  an  existing  method-of- moments  computer  code, 
thus  requiring  no  significant  code  modifications; 

.  accurately  constructs  the  surface  current,  impedance,  etc.; 

.  utilizes  simple  interpolation  functions  that  require  only  a  few  coefficients; 

.  apply  to  a  wide  variety  of  configurations; 
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.  impose  few  restrictions  on  the  geometry  modeling; 

.  are  independent  of  excitation. 

Each  element  of  the  matrix  is  approximated  by  the  ratio  of  two  polynomials  [1]  given  as 


CM)  = 


a  n  +  a  if  +  a  -,f2+...  +  a  fp 

mn,  0  mn,lJ  mn,2J  run ,  p  J 

+Kn,lf2  +-  +  Kn,dfd 


(4.1) 


where  /  denotes  the  frequency,  p  denotes  the  order  of  the  numerator  polynomial,  d 
denotes  the  order  of  the  denominator  polynomial,  and  amni  and  bmni  are  the  mn,h 
elements  of  the  complex  coefficient  matrices  [a,.]  and  [/?  ] ,  respectively.  When  the 
numerator  and  the  denominator  have  the  same  or  nearly  the  same  degree,  the  ratio  of 
polynomials  representation  of  matrix  element  Cmn  is  often  better  than  a  polynomial 
approximation  [10].  In  this  approach  p+d+1  coefficient  matrices  are  computed. 

The  Cmn  matrix  interpolation  begins  by  defining  the  order  of  the  polynomials  in 
(4.1).  Next,  p  +  d+1  frequencies  within  the  frequency  range  of  interest  are  chosen.  A 
N  xN  matrix  is  computed  at  each  selected  frequency  / .  The  elements  of  these  matrices 
are  then  used  to  determine  the  p  +  d+1  complex  coefficient  matrices  as  shown  in  (4.4). 

The  mn'1’  element  of  each  coefficient  matrix  is  then  substituted  into  (4.1)  to  compute  the 
corresponding  elements  of  the  Cmn  matrix  at  each  intermediate  frequency. 

The  Cmn  matrix  elements  are  highly  dependent  upon  the  resonant  characteristics 

of  the  structure.  The  choices  of  the  polynomial  orders  in  the  numerator  and  denominator 
depend  upon  the  proximity  of  adjacent  structure  resonances  and  the  overall  evaluation 
frequency  range.  Since  one  typically  does  not  know  the  location  of  the  resonances 
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beforehand,  there  is  no  convenient  way  to  quickly  determine  the  order  of  the  numerator 
and  denominator  polynomials.  In  this  work,  the  same  order  for  the  numerator  and  the 
denominator  polynomials  is  used.  Equation  (4.1)  can  be  rewritten  as 


CJf)+CJf)b„J  +  CJf)bm,i2f2+...  +  CJf)b„Jfd 
=  a  n  +a  if+a  7f2+...  +  a  fp.  (4.2) 

mn,  0  mn,±J  mn,2J  mn,pJ 

The  unknowns,  amni,  i  =  0,  ...,  p ,  and  bmni,  i  =  1,  ...,  d  ,  are  then  calculated  as  a  solution 
of  the  matrix  equation  where  /  is  the  selected  frequency  within  the  band  that  is  chosen 
for  direct  evaluation.  The  matrix  equation  is  of  the  form 

[a][x]  =  [b]  (4.3) 

where  [x]  is  the  unknown  column  matrix.  The  complete  equation  is  shown  in  (4.4). 


“1  fi  fi  fi  fi  -CmnJ i  -CmnJ2  -CmnJ?  -CmnJiff 

^mn ,  0 

Cmn.f i 

1  fi  fi  fi  ft  -cmn,hf2  -cmnJjl  -CmnJjl  -CmnJif* 

Cnm.fi 

1  . 

1  .  .  . 

^ mn,p-l 

If.. 

a 

mn,p 

= 

1  .  .  . 

b 

mn,L 

1  . 

!/«/«  ■  • 

i  n  P  fi  fi  -c, . -c„„,ji  -c, . ,j]  -cmfji 

b  . 

mn,a  J 

_  C mn .  fd  _ 

(4.4) 


The  flow  chart  in  Fig.  4-1  shows  the  polynomial  interpolation  method  as  applied 
to  a  narrow  slot  in  an  infinite  PEC  screen  backed  by  a  rectangular  cavity  with  a  thin-wire 
probe  in  it.  The  slot  and  the  probe  terms  in  the  presence  of  the  rectangular  cavity  are 


50 


computed  by  the  direct  M  oM  as  their  resonance  behavior  is  not  captured  with  a  10th- 
order  polynomial  interpolation  method  (see  Fig,  4-2)  and  it  is  advantageous  to  use  the 
direct  computation  instead.  The  slot  and  the  probe  terms  in  the  presence  of  the 
rectangular  cavity  vary  rapidly  with  frequency  and  hence  it  is  difficult  to  approximate 
these  terms  with  the  10th  order  interpolation  method.  The  coupling  terms  and  the  slot  self¬ 
terms  are,  however,  computed  by  the  use  of  the  ROPs  technique.  Eventually  all  the 
matrix  elements  are  assembled  together  for  impedance  calculations. 


Polynomial 

Interpolation 


Fig.  4-1.  M  atrix  interpolation  methodology  -  ROPs  technique  as  applied  to  a  narrow 
slot  in  an  infinite  PEC  screen  backed  by  a  rectangular  cavity 
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Fig.  4-2. 


Structure  Parameters 
2L  =12.5  mm 
2w  =  0.5  mm 
h  =  7.62  mm 
c  =47.275  mm 
d  =25.4  mm 


Frequency  (GHz) 

(a) 


(b) 


(a)  Slot  term  ( Y™ )  in  a  rectangular  cavity  for  the  structure  shown  in  Fig.  4-4. 

(b)  Probe  term  (Zn)  in  a  rectangular  cavity  for  the  structure  shown  in  Fig.  4- 
4.  (Single  matrix  element  computed  by  direct  M  oM  and  ROPs  interpolation) 
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3.  Results 

To  investigate  the  accuracy  of  the  matrix  interpolation  technique,  the  antenna 
input  impedance  is  illustrated  in  Fig.  4-3  for  a  dipole  antenna  of  length  L  =  0.5  m  over 
a  range  from  200  to  1400  MHz,  with  1st  -order,  2nd  -order  and  4th  -order  polynomial 
interpolations.  The  antenna  is  subdivided  into  20  equal  segments  in  length  of  about 
0.12A  at  1400  M  Hz  where  N  =  19  basis  functions  are  used.  The  results  show  close 
agreement  between  the  direct  MoM  and  2nd  and  4th-order  interpolations.  The  4th-order 
interpolation  curve  in  Fig.  4-3  provides  very  accurate  data  in  the  sense  that  the  error  in 
interpolating  [Zmn\  is  negligible.  However,  some  error  is  observable  in  the  2nd  and  re¬ 
order  cases.  The  deviation  from  the  direct  M  oM  results  around  640  M  Hz  and  1120  MHz 
is  a  good  measure  of  the  interpolation  error.  The  4th-order  interpolation  data  also  shows 
some  deviations,  but  they  are  much  smaller. 

One  may  use  the  polynomial  interpolation  technique  for  the  analysis  of  probes 
and  slots  in  rectangular  cavities,  the  results  of  which  are  discussed  in  detail  here.  As 
shown  in  Fig.  4-4,  the  structure  under  examination  is  a  narrow  slot  in  an  infinite  PEC 
screen  backed  by  a  rectangular  cavity  with  a  thin-wire  probe  present  inside  the  cavity. 
The  length  of  the  thin-wire  probe  is  h  and  its  radius  is  r.  The  probe  axis  is  located  in  the 
yz  plane  with  d  representing  its  displacement  from  the  z  =  0  plane.  The  slot  length  is 
2 L  and  its  width  is  2 w .  The  width,  height,  and  depth  of  the  cavity  shown  in  the  figure 
are  2 a,  2b  and  c,  respectively.  The  medium  inside  the  cavity  and  in  the  left  half-space  is 
characterized  by  {^e).  The  probe  is  connected  via  a  50Q  coaxial  cable  to  a  signal 
source  having  an  internal  impedance  of  50Q . 
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Fig,  4-5(a)  and  Fig,  4-5(b)  show  that  the  slot  self-coupling  elements  can  be 
interpolated  to  an  accurate  approximation  through  the  ratio  of  polynomials  interpolation 
method,  One  may  use  the  10th-order  ROPs  approximation  to  compute  the  slot  self¬ 
coupling  terms  and  Y2S3.  The  agreement  between  the  direct  M  oM  and  the 

interpolated  results  is  excel  lent  throughout  the  frequency  range  of  6  GHz  to  12  GHz. 
Fig.  4-6  and  Fig,  4-7  show  the  probe-to-slot  coupling  elements  G-  and  C- 

which  are  similar  to  C™  defined  in  (2,47).  One  again  sees  good  agreement  between  the 

direct  MoM  and  the  10th-order  polynomial  interpolated  results.  To  understand  and 
explain  these  results,  the  ROPs  representation  was  used  to  investigate  how  accurately  it 
can  compute  a  single  matrix  element.  The  coefficients  for  the  10th-order  ROPs 
approximation  were  computed  by  evaluating  matrix  elements  at  21  frequencies.  Fig.  4-6 
shows  that  6th-order  ROPs  capture  the  behavior  of  the  C^4  term  really  well  at  the  lower 

end  of  the  frequency  spectrum  but  does  not  at  the  higher  frequencies.  Fig.  4-8  shows  the 
si ot-to- probe  coupling  term  C*[,  which  is  similar  to  C™  defined  in  (2.52),  calculated  by 

direct  MoM  and  interpolation  techniques.  Again  the  agreement  is  excellent  throughout 
the  frequency  range. 

The  real  and  imaginary  parts  of  the  input  impedance  of  the  thin-wire  probe  in  the 
frequency  range  between  6  GHz  and  12  GHz  are  shown  in  Fig.  4-9  for  structure  A  of 
Table  4.1.  Numerical  results  are  obtained  by  using  the  direct  MoM  and  10th-order  ROPs 
interpolation.  Examining  the  figure,  the  TEm  mode  is  seen  to  be  excited  near  7.3  GHz 

in  both  methods,  which  is  very  close  to  the  resonant  frequency  of  this  mode  (7.5  GHz). 
The  TEm  mode  is  the  second  mode  excited  for  this  structure  with  its  excitation 
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occurring  near  10.3  GHz .  Since  the  probe  is  located  at  d  ^  c/2,  the  change  seen  in  the 
input  impedance  because  of  the  excitation  of  this  mode  is  not  as  drastic  as  the  excitation 
of  the  TE  mode.  The  direct  M  oM  results  and  interpolated  results  show  complete 

agreement  throughout  the  frequency  range. 

Fig.  4-10  shows  the  real  and  imaginary  parts  of  the  input  impedance  of  the  thin- 
wire  probe  for  structure  B  of  Table  4.1  in  the  frequency  range  between  6  GHz  and 
12  GHz.  Data  from  direct  MoM  is  compared  with  the  numerical  results  obtained  by 
using  10th-order  ROPs.  The  agreement  between  the  two  methods  is  very  good  over  the 
entire  frequency  range.  Examining  the  figure,  the  TEm  mode  is  excited  at  about 

7.27  GHz  in  both  the  methods  and  since  the  wire  is  located  at  d  ^  c/2 ,  the  TE[m  mode 
is  excited  slightly. 

Fig.  4-11  shows  the  cumulative  computation  time  for  direct  M  oM  and  that  for  the 
interpolation.  The  bulk  time  in  the  interpolated  method  is  utilized  in  the  calculation  of  the 
21  selected  frequencies  and  the  interpolation  procedure,  in  comparison,  is  a  fast  process. 
Examining  the  figure,  approximately  21000  seconds  were  required  for  the  computation 
of  matrix  elements  at  the  21  selected  frequencies  and  the  frequency  sweep  over  the  range 
from  6  GHz  to  12  GHz  was  a  much  faster  process.  In  contrast,  the  calculation  of  matrix 
elements  in  the  direct  method  is  a  computationally  very  heavy  task.  The  time  advantage 
in  the  ROPs  method  is  apparent  from  the  figure. 

The  [Y]  matrix  interpolation  and  coupling  matrix  interpolation  methodologies 
have  been  used  with  the  method  of  moments  in  order  to  significantly  reduce  the 
computation  time  required  for  the  wideband  performance  evaluation  of  slot-coupled 
cascaded  rectangular  cavities.  These  interpolation  methods  do  not  require  a  prior 
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knowledge  of  cavity  resonance  characteristics.  The  results  included  a  comparison  of  the 
time  savings  of  the  interpolation  method  over  the  direct  MoM  approach.  Sometimes 
large-order  polynomials  are  required  to  accurately  reproduce  the  elements  via  matrix 
interpolation  and  still  they  have  a  good  time  advantage  over  the  direct  method.  Higher- 
order  polynomials  for  matrix  interpolation  require  the  storage  of  more  coefficients  than 
the  lower-order  matrix  interpolation  techniques  and  require  more  computation  time.  The 
advantage  of  a  particular  interpolation  method  such  as  ROPs  technique  to  simple 
interpolation  techniques  such  as  the  quadratic  interpolation  is  based  upon  the  general 
resonant  nature  of  the  cavity  to  be  modeled.  When  many  resonances  occur  within  a 
particular  frequency  range,  ROPs  matrix  interpolation  is  faster  and  requires  less 
computational  resources  than  the  direct  method.  The  results  presented  in  this  chapter 
demonstrate  the  applicability  of  the  ROPs  matrix  interpolation  methods. 
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Antenna  Parameters 
Length  of  dipole  =  0.5  m 
Radius  of  dipole  =  0.0025  m 


(a) 


(b) 

Fig.  4-3.  Dipole  antenna  input  impedance  using  polynomial  interpolation  technique, 
(a)  real  part  of  impedance,  (b)  imaginary  part  of  impedance. 
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Table  4.1 

Structure  Parameters 


Parameter 

Structure  A 

Structure  B 

Cavity  width,  2 a 

22.86 

22.86 

Cavity  height,  2 b 

10.16 

10.16 

Cavity  depth,  c 

47.625 

41.275 

Slot  length,  2 L 

12.5 

15.0 

Slot  width,  2w 

0.5 

0.5 

Probe  length,  h 

7.62 

6.35 

Probe  radius,  r 

0.254 

0.254 

Probe  position,  d 

25.4 

17.4625 

*  A 1 1  parameters  are  in  millimeters 


Table  4.2 

Resonant  Frequencies" 


Structure  A 

Structure  B 

i— i 

o 

i— i i 

7.27 

7.50 

TE102 

9.09 

9.79 

TEm 

11.50 

12.72 

TE201 

13.49 

13.61 

TEm 

14,20 

15.94 

Probe  half-wave  resonance 

9.84 

11.81 

Slot  half-wave  resonance 

12.00 

10.00 

*  A 1 1  frequencies  are  in  GHz. 


Real  part  of  S(1,1 )  &  S(2,3) 
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Structure  A  Parameters 

yt 

2L  =12.5  mm 

PEC  walls 

‘ 

2w  =  0.5  mm 

h  =  7.62  mm 

r . d  "Y 

c  =47.275  mm 

/  * 

d  =25.4  mm 

P 

Probe  \  \ 

\\ 

-C  \\  Z 

Screen 


(a)  (b) 

Fig.  4-5.  Slot  self-coupling  elements-  Y°v  Y, 's3  -  (a)  real  part  (b)  imaginary  part 


Frequency  (GFIz) 

Fig,  4-6.  Coupling  element  from  probe  to  slot  -  C*p4  (single  matrix  element  computed 
by  direct  M  oM  and  ROPs  interpolation) 
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Fig.  4-7  Coupling  element  from  probe  to  slot  -  CS2P6  (single  matrix  element  computed 
by  direct  M  oM  and  ROPs  interpolation) 


Frequency  (GFIz) 

Fig.  4-8  Coupling  element  from  slot  to  probe  -  CPA\  (single  matrix  element  computed 
by  direct  M  oM  and  ROPs  interpolation) 
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Structure  A  Parameters 
2L  =12.5  mm 
2w  =  0.5  mm 
h  =  7.62  mm 
c  =  47.625  mm 
d  =25.4  mm 


Frequency  (GHz) 

(a) 


(b) 


Fig.  4-9.  Input  impedance  of  the  thin-wire  probe  inside  a  slotted  cavity,  (a)  real  part  of 
input  impedance,  (b)  imaginary  part  of  input  impedance.  (30  frequency  points 
computed  for  direct  M  oM  ) 
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Structure  B  Parameters 
2L  =15.0  mm 
2w  =  0.5  mm 
h  =  6.35  mm 
c  =41.275  mm 
d  =  17.4625  mm 


Frequency  (GHz) 

(a) 


(b) 


Fig.  4-10.  Input  impedance  of  the  thin-wire  probe  inside  a  slotted  cavity,  (a)  real  part  of 
input  impedance,  (b)  imaginary  part  of  input  impedance.  (30  frequency  points 
computed  for  direct  M  oM  ) 
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★  Direct  MoM  -  Cpu  time 
-  10th  order  Interpolation  -  Cpu  time 


Frequency(GHz) 


Fig,  4-11.  Cumulative  CPU  time  for  direct  method  (MoM)  and  polynomial 
interpolation  method  for  structure  A.  (60  frequency  points  computed  for 
direct  MoM  and  ROPs  interpolation) 
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4.  Conclusions 

In  this  study  an  investigation  of  an  efficient  method  for  the  computation  of  the 
broadband  characteristics  of  slot-coupled  cascaded  rectangular  cavities  using  frequency 
interpolation  of  MoM  matrix  is  presented.  The  ROPs  method  of  MoM  matrix 
interpolation  significantly  reduced  the  computation  time  required  for  the  wideband 
performance  evaluation  of  cavity  responses.  Good  agreement  showing  the  efficiency  and 
the  accuracy  of  the  ROPs  method  is  observed  in  the  results  obtained  by  the  ROPs  matrix 
interpolation  method  and  the  results  obtained  by  the  direct  M  oM  .  Data  included  the  CPU 
time  comparison  between  the  two  methods. 

Appendix:  Use  Of  X eon  M  ulti-Processorsto  Reduce  CPU  Time 
The  numerical  solution  of  the  cavity  problem  involves  computation  of  very 
slowly  convergent  infinite  series  in  the  calculation  of  the  matrix  elements  occurring  in  the 
method  of  moments  (MoM)  solution  of  integral  equations,  hence  the  enormous  time 
required  to  compute  the  matrix  elements.  This  led  to  the  farming  of  codes  on  Xeon  multi¬ 
processors,  wherein  the  original  FORTRAN  code  is  fed  into  different  Xeon  processors 
simultaneously  through  the  use  of  PBS  (portable  batch  system)  scripts.  The  main  script 
controls  the  feeding  of  the  sub-codes  through  the  queue  command  and  also  assembles  the 
results  at  the  end  of  the  runs.  The  full  process  is  illustrated  in  the  flow  chart  in  Fig.  4-12. 
The  computationally  heavy  segment  of  the  FORTRAN  code  involves  the  calculation  of 
matrix  elements.  The  main  script  submits  the  '  n '  subordinate  scripts  one  after  the  other 
into  the  queue  and  the  subordinate  scripts  are  submitted  to  the  different  nodes  for  separate 
calculations.  The  main  script  also  assembles  all  the  resultant  elements  into  a  complete 
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matrix  file.  An  important  criterion  is  the  synchronization  of  all  '  n  '  sub-scripts  such  that 
the  next  stage  of  computation  starts  after  all  '  n  '  codes  exit  successfully.  The  entire  (n,n) 
matrix  file  is  depicted  in  Fig.  4-13  along  with  the  method  of  splitting  them  into  21 
(«  =  21)  different  segments  to  be  fed  into  various  nodes.  The  Yx  and  Y{  terms  in  the 
figure  are  the  admittance  terms  due  to  slot  1,  and,  Y2  and  Y{  are  the  admittance  terms 
due  to  slot  2.  The  elements  of  the  sub-matrices  representing  the  coupling  from  the  probe 
to  slot  1  are  given  by  Cslp  and  the  coupling  from  probe  to  slot  2  is  given  by  Cs2p .  Cpsl 
represents  the  slot  1  to  probe  coupling  term  and  Cpsl  represents  the  slot  2  to  probe 
coupling  term,  Z  is  the  impedance  sub-matrix  term  for  the  probe. 

The  CPU  time  results,  Fig.  4-14,  for  the  performance  of  parallel  codes  in  Xeon 
multi-processors  are  compared  with  a  run  in  an  Alpha  system  of  dual  processors  (serial) 
for  structure  A  parameters.  The  results  show  a  massive  improvement  in  cumulative  CPU 
time  when  run  on  Xeons  due  to  the  farming  across  multiprocessors,  an  improvement  of 
nearly  94%  from  the  serial  case  for  a  wideband  evaluation  between  6  GHz  and  12  GHz . 
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Fig.  4-12.  Farming  of  FORTRAN  program  across X eon  multi-processors 
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Structure  Geometry 


Fig.  4-13.  M  atrix  fill  mechanism  for  a  single  rectangular  cavity  with  two  slots  and  one 
probe. 


CD 

C/D 

LO 

o 


CD 

E 

\— 

ZD 

Q_ 

o 

0 

> 

To 

E 

=3 

o 


Alpha  Time  (Serial)  Vs 

Xeon  Time  (Split  Parallel) 

X  Split  Parallel 
o  Serial 

Q 

- 

O 

O 

O 

O 

O 

O 

O 

Q 

Q 

O 

Q 

O 

O 

o 

o 

o 

-  ,  v ,  v  ,vi  x  ,x,  X  |X,  X 

X 

X 

X 

X 

X 

X 

x” 

X 

o  x 

o  x- 

6  7  8  9  10  11  12 


Frequency  (GHz) 


Fig.  4-14.  Comparison  of  cumulative  CPU  time  for  Alphas  and  Xeon  processors 
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Chapter  5 


Numerical  Analysis  Of  Slot-Coupled  Cascaded 

Rectangular  Cavities 

Numerical  results  for  signal  propagation  through  cascaded  rectangular  cavities  are 
presented  in  this  chapter.  The  cascaded  rectangular  cavities  are  interconnected  by  narrow 
slots  and  have  thin-wire  probes/posts  present  in  them.  The  input  impedance  at  the  center 
of  the  slot,  input  impedance  at  the  base  of  the  thin-wire  probe  and  the  electric  field 
shielding  (EFS)  effectiveness  [6]  inside  the  rectangular  cavity  are  calculated  for  single, 
double,  and  triply  cascaded  rectangular  cavities. 

1.  Introduction 

A  numerical  analysis  of  thin-wire  probes  inside  slotted  cascaded  rectangular 
cavities  is  undertaken.  The  method  of  analysis  is  based  on  solving  the  coupled  integral 
equations  developed  for  the  electric  current  on  the  probes  and  the  equivalent  magnetic 
current  in  the  slots  (see  Chapter  2).  The  structure  considered  for  investigation  is  a  narrow 
slot  in  an  infinite  PEC  screen  backed  by  cascaded  rectangular  cavities  that  contain  thin- 
wire  probes/posts  in  them  and  which  are  interconnected  by  narrow  slots.  The  analysis  of 
this  structure  requires  coupled  integral  equations  which  include  the  free-space  periodic 
Green’s  function  as  part  of  their  kernels.  The  coupled  integral  equations  are  solved  by  the 


method  of  moments  [2]. 
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The  EFS  [6]  is  calculated  at  the  center  of  each  rectangular  cavity  to  evaluate  the 
effect  of  the  transmission  path/environment  on  the  entering  transient  signal  and  is  defined 
as 


EFS  (dB)  =  -20  log 


Ey(x,y,z ) 
EyC(x,y,z) 


(5.1) 


where  E  is  the  amplitude  of  the  electric  field  at  a  specified  point  in  the  presence  of  the 
enclosure  and  E‘"c  is  the  amplitude  of  the  electric  field  at  the  same  point  with  the 
enclosures  removed.  The  field  E"u:  is  given  by 


=  V  e~jkz 


(5.2) 


where  V0  is  the  input  voltage  at  the  center  of  the  slot  and  E  is  obtained  from 


Ey(x,  y,z)  =  — 
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(5.3) 


where  m(x)  is  the  axial  variation  of  the  magnetic  current  on  the  slot  and  I(y)  is  the 
electric  current  on  the  thin-wire  probe.  The  term  K(y,y')  is  similar  to  (y,yr)  in 

(2.35)  and  Gd(y,x',y ')  is  similar  to  Gd  {y,x' ,y')  in  (2.25).  Equation  (5.3)  is  applicable 
for  a  rectangular  cavity  with  a  single  slot  and  probe.  If  the  rectangular  cavity  has  two 


slots,  then 
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Ey(x,  y,z)  =  — 
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where  m,  (x)  and  m2(x)  are  the  axial  variations  of  the  magnetic  currents  on  slot  1  and 
slot  2  and  K'‘v(y,  y') ,  GJ1  (y,x',  y')  and  Gd2(y,x',y')  are  similar  to  terms  defined  in 
(2.19). 

The  expression  for  the  electric  field  on  the  non-illuminated  side  (Fig.  5-1)  of  the 
final  cavity  screen  is  given  below.  The  magnetic  current  in  the  last  slot  is  determined  by 
the  use  of  the  MoM  solution  procedure  and  one  can  compute  the  electric  field,  E  , 

behind  the  final  cavity  screen  from 


p  ,  ,17,  r/f  1/^r 

Ey(W)=  —  m-(X)  I  K 

Ln  -L,  y'=-wy]W  -y 


where  ‘  n  ’  is  the  index  of  the  slot  current  on  the  end  wall.  By  expanding  mn  (x)  in  pulse 


basis  functions,  one  can  write  E  as 
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Fig.  5-1.  Two  cascaded  rectangular  cavities  showing  the  non-illuminated  point  behind 
the  last  cavity  screen. 


As  an  investigation  of  the  method  devised  to  calculate  EFS  of  a  shielded  enclosure,  the 
shielding  effectiveness  of  a  single  rectangular  cavity  with  a  narrow  slot  excited  by  a 
plane-wave  source  was  computed  and  EFS  plotted  at  the  center  of  the  cavity  for  a 
frequency  range  of  100  to  1000  MHz .  The  EFS  plot  was  checked  against  results  in  [4] 
and  [6]  and  found  to  agree  well.  A  circuital  approach  to  predict  the  field  distributions  on 
rectangular  apertures  backed  by  rectangular  cavities  is  discussed  by  Azaro  in  [4]  and  the 
multilevel  fast  multipole  method  (MLFMM)  algorithm  is  used  for  the  computation  of  the 


EFS  parameters  by  Volakis  in  [6]. 
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(■ H‘,E‘ ) 

Slot 


PEC 


(b) 


Fig.  5-2.  Shielding  effectiveness  at  the  center  of  the  rectangular  cavity.  Enclosure 
dimensions:  30  x  12  x  30  (cm),  (a)  Aperture  dimensions:  20  x  3.0  (cm) 
(b)  Aperture  dimensions:  10  x  0.5  (cm) 
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The  excitation  used  to  illuminate  the  slot  is  shown  in  Fig.  5-3(a).  It  is  a  differentiated 
Gaussian  pulse  and  its  Fourier  transform  is  shown  in  Fig.  5-3(b).  The  energy  is  centered 
around  1  GHz  and  only  the  positive  frequencies  are  shown  in  the  figure. 


(a) 


Fig.  5-3. 


(a)  Differentiated  Gaussian  pulse  excitation:  (n  —  n 0)exp( - - 

It, 

t  =  67.45  secs  and  no  =  3r  .  (b)  FFT  spectrum  (magnitude)  of  the 
differentiated  Gaussian  pulse  with  its  energy  centered  around  1  GHz. 


)  with 
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2.  Single  Rectangular  Cavity  With  A  Thin- Wire  Probe 
And  A  Narrow  Slot 

As  shown  in  Fig.  5 -4(a),  the  structure  under  examination  is  a  narrow  slot  in  an 
infinite  PEC  screen  backed  by  a  rectangular  cavity  with  a  thin-wire  probe/post  in  it.  The 
side  view  of  the  structure  is  illustrated  in  Fig.  5-4(b).  As  seen  from  the  figure,  it  is  a 
single  rectangular  cavity  where  the  length  of  the  thin-wire  probe  is  h  with  its  radius 
specified  as  r .  The  probe  axis  is  located  in  the  yz  plane  with  d  representing  its 
displacement  from  the  z  =  0  plane.  The  slot  length  is  2 L  and  the  width  is  2  w .  The 
width,  height,  and  depth  of  the  cavity  shown  in  the  figure  is  2 a ,  2b  and  C  ,  respectively. 
The  medium  inside  and  in  the  left  half-space  is  characterized  by  ( ju ,  s) . 

The  numerical  results  in  the  following  pages  show  the  input  impedance  of  the  slot 
and  the  probe  for  different  structures.  The  input  impedance  is  calculated  from  the 
magnetic  currents  and  the  electric  currents  which  are  obtained  by  solving  the  matrix 
equation  as  given  in  (2.37).  The  orientation  of  the  thin-wire  probe  in  the  cavity  and  the 
orientation  of  the  slot  in  the  cavity  do  not  allow  the  TE20p  and  TE01p  modes  to  be  excited 

in  the  cavity.  One  can  find  the  modes  excited  by  calculating  the  resonant  frequency  or  the 
wavelength  of  each  mode  and  examining  the  currents  [5].  One  finds  that  the  TEm  mode 

is  the  lo west-order  mode  excited  for  each  structure. 

The  real  and  imaginary  parts  of  the  input  impedance  at  the  center  of  the  slot  for 
Structure  A  are  shown  in  Fig.  5-5.  The  upper  limit  of  the  frequency  spectrum  is  5  GHz . 
In  the  figure,  one  can  see  that  the  lowest-order  TEm  mode  is  excited  at  about  3.78  GHz 

which  is  very  close  to  the  resonant  frequency  of  this  mode  (3.40  GHz),  as  shown  in 
Table  5.2.  The  resonance  associated  with  the  thin-wire  probe  is  seen  at  about  2.58  GHz  , 
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which  is  close  to  the  resonance  associated  with  the  quarter-wave  monopole  (3.125  GHz ), 
with  a  zero  crossing  in  the  imaginary  part  of  the  input  impedance.  The  narrow-slot 
resonance  is  seen  at  about  4.19  GHz  which  is  very  close  to  the  calculated  resonance  of 
4.28  GHz . 

Fig.  5-6  shows  the  real  and  imaginary  parts  of  the  input  impedance  computed  at 
the  center  of  the  slot  for  Structure  B  with  an  upper  frequency  limit  of  5  GHz  .  Examining 
the  figure  reveals  that  the  lowest-order  TEm  mode  is  excited  at  about  3.27  GHz  which  is 

very  close  to  the  resonant  frequency  of  this  mode  (3.40  GHz )  as  shown  in  Table  5.2.  The 
figure  also  shows  that  the  resonance  associated  with  the  narrow  slot  is  seen  at  about 
4.37  GHz  which  is  close  to  the  calculated  resonance  of  4.28  GHz  . 

Table  5.3  shows  parameters  for  Structure  C  for  which  the  input  impedance  at  the 
center  of  the  slot  and  at  the  base  of  the  thin-wire  probe  is  computed  by  driving  them 
separately.  The  slot  is  excited  with  the  frequency  spectrum  obtained  from  the  FFT  of  a 
differentiated  Gaussian  pulse  as  described  in  Fig.  5-3(b)  and  the  probe  is  fed  with  a  1.0  V 
source  at  its  base.  The  EFS  for  the  structure  is  computed  at  the  center  of  the  cavity  using 
(5.1).  In  the  Fig.  5-7,  one  can  see  that  the  lowest-order  mode,  TEm,  is  excited  at  about 

3.41  GHz,  which  is  very  close  to  the  resonant  frequency  of  this  mode  (3.40  GHz),  as 
shown  in  Table  5.4.  The  resonance  associated  with  the  thin-wire  probe  is  seen  at  about 
2.44  GHz ,  which  is  close  to  the  resonance  associated  with  the  quarter-wave  monopole 
( 2.67  GHz ),  with  a  zero  crossing  in  the  imaginary  part  of  the  input  impedance.  The 
narrow-slot  resonance  is  seen  at  about  4.29  GHz  which  is  extremely  close  to  the 
calculated  resonance  of  4.28  GHz .  The  input  impedance  at  the  base  of  the  probe  is 
shown  in  Fig.  5-8.  The  resonances  associated  with  the  cavity,  slot  and  the  probe  are  very 
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clearly  defined  at  the  same  frequencies  as  seen  for  the  slot  impedance.  The  EFS  is 
computed  at  the  center  of  the  cavity  using  (5.1).  The  plot  (Fig.  5-9)  shows  that  EFS 
decreases  in  the  neighborhood  of  the  resonances  as  expected. 

Fig.  5-10  shows  the  real  and  imaginary  parts  of  the  input  impedance  computed  at 
the  center  of  the  slot  for  Structure  D  with  an  upper  frequency  limit  of  5  GHz  .  The  same 
driving  mechanism  is  used  as  described  in  the  previous  section.  As  the  figure  shows,  the 
lowest-order  mode,  TEm,  is  excited  at  about  4.22  GHz  which  is  close  to  the  resonant 
frequency  of  this  mode  (3.86  GHz),  as  shown  in  Table  5.4.  The  figure  also  shows  that 
the  resonance  associated  with  the  narrow  slot  is  seen  at  about  3.01  GHz  which  is 
extremely  close  to  the  calculated  resonance  of  3.0  GHz .  The  resonance  associated  with 
the  thin- wire  probe  is  seen  at  about  1.84  GHz ,  which  is  close  to  the  resonance  associated 
with  the  quarter-wave  monopole  (2.08  GHz).  The  EFS  plot  (Fig.  5-11)  shows  that  the 
shielding  effectiveness  decreases  near  1.84  GHz,  3.11GHz,  and  4.24  GHz  which 
corresponds  to  the  three  resonances  captured  in  the  impedance  plots. 
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Table  5.1 

Structure  Parameters 


Parameter 

Structure  A 

Structure  B 

Cavity  width,  2 a 

60.0 

60.0 

Cavity  height,  2b 

34.0 

32.0 

Cavity  depth,  c 

65.0 

65.0 

Slot  length,  2 L 

35.0 

35.0 

Slot  width,  2 w 

0.5 

0.5 

Probe  length,  h 

24.0 

12.0 

Probe  radius,  r 

0.25 

0.25 

Probe  position,  d 

20.0 

20.0 

*A11  parameters  are  in  millimeters 


Table  5.2 

Resonant  Frequencies" 


Structure  A 

Structure  B 

TEm 

3.40 

3.40 

TE\()2 

5.25 

5.25 

TE103 

7.36 

7.36 

TE20l 

5.5 

5.5 

9.56 

9.56 

Probe  half-wave  resonance 

3.125 

6.25 

Slot  half-wave  resonance 

4.28 

4.28 

*A11  frequencies  are  in  GHz. 
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Fig.  5-4.  A  narrow  slot  in  an  infinite 
thin- wire  probe  in  it.  (a)  pet 
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Fig.  5-5.  Input  impedance  at  the  center  of  the  slot,  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Fig.  5-6.  Input  impedance  at  the  center  of  the  slot,  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Table  5.3 

Structure  Parameters 


Parameter 

Structure  C  (mm) 

Structure  D  (mm) 

Cavity  width,  2 a 

60.0 

53.0 

Cavity  height,  2b 

34.0 

40.0 

Cavity  depth,  c 

65.0 

58.0 

Slot  length,  2  L 

35.0 

50.0 

Slot  width,  2 w 

0.5 

0.5 

Probe  length,  h 

28.0 

36.0 

Probe  radius,  r 

0.25 

0.25 

Probe  position,  d 

4.0 

20.0 

Table  5.4 

Resonant  Frequencies 


Structure  C 

Structure  D 

T^iot 

3.40 

3.86 

teXQ2 

5.25 

5.87 

TEm 

7.36 

7.36 

EE  201 

5.5 

5.5 

EE  104 

9.56 

9.56 

Probe  half-wave  resonance 

2.67 

2.08 

Slot  half-wave  resonance 

4.28 

3.00 

**  All  Frequencies  are  in  GHz 
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Structure  C  Parameters 
2L  =35.0  mm 
2w  =  0.5  mm 
h  =  28.0  mm 
c  =  65.0  mm 
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Fig.  5-7.  Input  impedance  at  the  center  of  the  slot,  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Structure  C  Parameters 
2L  =35.0  mm 
2w  =  0.5  mm 
h  =  28.0  mm 
c  =  65.0  mm 
d  =  4.0  mm 


(a) 


Fig.  5-8.  Input  impedance  at  the  base  of  the  probe,  (a)  Real  part  of  impedance  (b) 
Imaginary  part  of  impedance. 
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Fig.  5-9.  EFS  at  the  center  of  the  rectangular  cavity  for  Structure  C. 
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Fig.  5-10.  Input  impedance  at  the  center  of  the  slot,  (a)  Real  part  of  impedance 
(b)  Imaginary  part  of  impedance. 


87 


Slot 

<Ji,s) 


(M>e) 


c/2 


PEC 


Fig.  5-11.  EFS  at  the  center  of  the  rectangular  cavity  for  Structure  D. 
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3.  Two  Cascaded  Rectangular  Cavities  Interconnected 
By  A  Narrow  Slot  With  Thin-Wire  Probes  Inside 

As  shown  in  Fig.  5-12,  the  structure  under  examination  is  a  narrow  slot  in  an 
infinite  PEC  screen  backed  by  two  cascaded  rectangular  cavities  each  having  a  thin-wire 
probe/post  inside  and  which  are  interconnected  by  a  narrow  slot.  The  lengths  of  the  thin- 
wire  probes  are  /?,  and  h2  with  their  radii  specified  as  r .  The  probe  axes  are  located  in 
the  yz  plane  with  d2  and  c  +  d2  representing  their  displacement  from  the  z  =  0  plane. 
The  slot  lengths  are  2Lj ,  2 L2  and  2 L3  respectively  and  their  widths  are  2 w  for  all  the 
slots.  The  width,  height,  and  depth  of  each  cavity  shown  in  the  figure  is  2 a  ,  2b  and  c . 
The  medium  inside  each  cavity,  in  the  left  half-space  where  z  <  0  and  in  the  right-half 
space  where  z>2c  is  characterized  by  (ju,s) . 

The  numerical  procedure  was  first  tested  by  computing  the  electric  and  magnetic 
currents  on  the  probes  and  the  slots  for  a  symmetrical  structure  which  is  excited  at  the 
second  (middle)  slot.  The  currents  are  shown  in  Fig.  5-13.  The  magnetic  current  on  slot  1 
and  slot  3  and  the  electric  current  on  probe  1  and  probe  2  are  plotted.  From  the  figure  it  is 
seen  that  the  agreement  between  the  slot  currents  and  the  probe  currents  is  very  good  at 
the  given  frequency,  which  is  to  be  expected.  Though  this  method  is  not  a  foolproof 
validation  of  the  numerical  process,  it  serves  as  a  partial  validation  of  the  numerical 
results. 

Table  5.4  shows  parameters  for  Structure  E  for  which  input  impedance  at  the 
center  of  slot  1  and  at  the  base  of  probe  1  in  the  first  cavity  is  computed  by  driving  them 
separately.  The  frequency  spectrum  obtained  from  the  FFT  of  a  differentiated  Gaussian 
pulse  as  defined  in  Fig.  5-3(b)  is  used  to  excite  the  narrow  slot  and  the  probe  is  fed  with  a 
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1.0  V  source  at  its  base.  The  EFS  for  the  structure  is  computed  at  the  center  of  each 
cavity  using  (5.1).  In  the  Fig.  5-14,  one  can  see  that  the  lowest-order  mode,  TEm,  is 
excited  at  about  3.38  GHz  which  is  extremely  close  to  the  resonant  frequency  of  this 
mode  (3.40  GHz ),  as  shown  in  Table  5.4.  The  resonance  associated  with  the  thin-wire 
probe  is  seen  at  about  2.43  GHz ,  which  is  close  to  the  resonance  associated  with  the 
quarter-wave  monopole  ( 2.67  GHz  ),  where  a  zero  crossing  in  the  imaginary  part  of  the 
input  impedance  is  seen.  The  narrow-slot  resonance  is  seen  at  about  4.26  GHz  ,  which  is 
extremely  close  to  the  calculated  resonance  of  4.28  GHz  .  The  resonance  of  the  second 
slot  (slot  2)  is  seen  at  4.41  GHz  and  the  resonance  of  the  third  slot  (slot  3),  though  not 
properly  excited,  exists  at  around  4.58  GHz .  The  input  impedance  at  the  base  of  the 
probe  is  shown  in  Fig.  5-15.  The  resonances  associated  with  the  cavity,  slot,  and  the 
probe  are  very  clearly  defined  at  the  same  frequencies  as  seen  for  the  slot  impedance.  The 
EFS  is  computed  at  the  center  of  each  cavity  (z  =  c/ 2,  z  =  3c/2  and  z  =  5c/2 )  using 
(5.1).  Equation  (5.6)  is  used  to  compute  E  at  z  =  5c/2  .  The  plot  shows  that  EFS  (Fig. 

5-16)  decreases  in  the  neighborhood  of  the  resonances,  as  expected. 

Fig.  5-17  shows  the  real  and  imaginary  parts  of  the  input  impedance  computed  at 
the  center  of  the  slot  for  Structure  F  with  an  upper  frequency  limit  of  5  GHz  .  Examining 
the  figure  reveals  that  the  lowest-order  mode,  TEm,  is  excited  at  about  4.22  GHz  which 
is  close  to  the  resonant  frequency  of  this  mode  (3.86  GHz),  as  shown  in  Table  5.4.  The 
resonance  of  the  second  cavity  is  seen  to  be  excited  slightly  at  about  4.35  GHz .  The 
figure  also  shows  that  the  resonance  associated  with  the  narrow  slot  is  seen  at  3.01  GHz , 
which  is  extremely  close  to  the  calculated  resonance  of  3.0  GHz  .  The  resonance  of  slot  2 
is  seen  at  3.11  GHz  .  The  resonance  of  slot  3  is  not  properly  excited  and  hence  is  not  seen 
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in  the  figure.  The  resonance  associated  with  the  thin- wire  probe  is  seen  at  about 
1.84  GHz,  which  is  close  to  the  resonance  associated  with  the  quarter-wave  monopole 
(2.08  GHz).  The  EFS  plot  in  Fig.  5-18  shows  that  the  shielding  effectiveness  decreases 
near  1.84  GHz  ,  3.11  GHz  ,  and  4.22  GHz  ,  which  corresponds  to  the  three  resonances 
captured  in  the  impedance  plots. 
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Table  5.3 

Structure  Parameters 


Parameter 

Structure  E  (mm) 

Structure  F  (mm) 

Cavity  width,  2 a 

60.0 

53.0 

Cavity  height,  2b 

34.0 

40.0 

Cavity  depth,  c 

65.0 

58.0 

Slot  length,  2LV  2 L2,  2 L, 

35.0 

50.0 

Slot  width,  2 w 

0.5 

0.5 

Probe  length,  1\ ,  h2 

28.0 

36.0 

Probe  radius,  r 

0.25 

0.25 

Probe  position,  dx ,  d2 

4.0 

20.0 

Table  5.4 

Resonant  Frequencies 


Structure  E 

Structure  F 

77?10i 

3.40 

3.86 

77?  102 

5.25 

5.87 

77?  103 

7.36 

7.36 

7^201 

5.5 

5.5 

7^104 

9.56 

9.56 

Probe  half-wave  resonance 

2.67 

2.08 

Slot  half-wave  resonance 

4.28 

3.00 

**  All  Frequencies  are  in  GHz 
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z  =  0  z  =  2c 
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Fig.  5-12.  Two  cascaded  rectangular  cavities  interconnected  by  a  narrow  slot  and  with 
thin- wire  probes  in  them,  (a)  perspective  view,  (b)  side  view. 
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Structure  Parameters 

2Lj  =  2  L7  =  2L3=35.0mm 

2w=  0.5  mm 

hl  =  /?,  =  24.0  mm 

C  =  65.0  mm 

2 a  =  60.0  mm,  2b  =  34.0  mm 
dl  =  d2  =  c/2 
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Fig.  5-13.  Magnetic  currents  in  the  slots  and  electric  currents  on  the  probes  due  to  a  1.0 
V  excitation  on  slot  2  (a)  Magnetic  currents  in  slot  1  and  slot  3  (b)  Electric 
currents  on  probe  1  and  probe  2. 


Structure  E  Parameters 

2L,  =  2 =  2 L3  =  35.0  mm 

2w  =  0.5  mm 

l\  =  =  28.0  mm 

c  =  65.0  mm 

cl{  =  =  4.0  mm 
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Fig.  5-14.  Input  impedance  at  the  center  of  slot  1.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Structure  E  Parameters 

2L,  =  2 Ln  =  2 L3  =  35.0  mm 

2 w=  0.5  mm 

h{  =  =  28.0  mm 

c  =  65.0  mm 

dt  =  d2  =  4.0  mm 


y 


Fig.  5-15.  Input  impedance  at  the  base  of  probe  1.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Fig.  5-16.  EFS  for  Structure  E. 


97 


Structure  F  Parameters 

2Lj  =  2 L,  =  2L3  =  50.0  mm 

2w=  0.5  mm 

/tj  =  /t0  =  36.0  mm 

c  =  58.0  mm 

<7,  =  <7,  =  20.0  mm 


-TvAA->  Slot  1 


(/".«)  .b 


<ju,e) 


(p,e) 


r 


2c 
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Fig.  5-17.  Input  impedance  at  the  center  of  slot  1.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Fig.  5-18.  EFS  for  Structure  F 
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4.  Three  Cascaded  Rectangular  Cavities  Interconnected 
By  Narrow  Slots  With  Thin- Wire  Probes  Inside 

As  shown  in  Fig.  5-19,  the  structure  under  examination  is  a  narrow  slot  in  an 
infinite  PEC  screen  backed  by  triply  cascaded  rectangular  cavities,  each  with  a  thin-wire 
probe/post  in  it  and  which  are  interconnected  by  narrow  slots.  As  seen  from  the  figure, 
the  structure  is  a  triply  cascaded  rectangular  cavity  where  the  lengths  of  the  thin-wire 
probes  are  \ ,  h2  and  h{  with  their  radii  specified  as  r .  The  probe  axes  are  located  in  the 

yz  plane  with  dl ,  c  +  d2  and  2c  +  d3  representing  their  displacement  from  the  z  =  0 
plane.  The  slot  lengths  are  2 L, ,  2 L2 ,  2 L3  and  2L4  respectively  and  their  widths  are  2w 
for  all  the  slots.  The  width,  height,  and  depth  of  each  of  the  cavity  shown  in  the  figure  is 
2a ,  2b  and  c ,  respectively.  The  medium  inside  each  cavity,  in  the  left  half-space  and  in 
the  right-half  space,  is  characterized  by  (//,  s) . 

The  numerical  procedure  was  first  tested  by  computing  the  electric  and  magnetic 
currents  on  the  probes  and  the  slots  for  a  symmetrical  structure  which  was  excited  at  the 
second  (middle)  probe.  The  currents  are  shown  in  Fig.  5-20.  The  equivalent  magnetic 
currents  in  slot  1  and  slot  4  are  plotted  in  Fig.  5-20(a)  and  in  slot  2  and  slot  3  are  plotted 
in  Fig.  5-20(b).  The  electric  currents  on  probe  1  and  probe  3  are  plotted  in  Fig.  5-20(c). 
From  the  figure  it  is  seen  that  the  agreement  between  the  slot  currents  and  the  probe 
currents  is  very  good  at  the  given  frequency  which  is  to  be  expected.  As  stated  earlier, 
this  method  enables  us  to  partially  validate  the  numerical  results  even  if  not  completely 
endorsing  it. 

Table  5.5  shows  parameters  for  Structure  G  for  which  input  impedance  at  the 
center  of  slot  1  and  at  the  base  of  probe  2  in  the  second  cavity  is  computed  by  driving 
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them  separately.  The  slot  is  excited  by  the  frequency  spectrum  obtained  from  the  FFT  of 
the  differentiated  Gaussian  pulse  (Fig.  5-3(b))  and  the  probe  is  driven  by  a  1.0  V  source 
applied  at  its  base.  The  EFS  for  the  structure  is  computed  at  the  center  of  each  cavity  to 
predict  the  influence  of  the  transmission  path  on  the  transient  signal.  One  can  see  from 
the  Fig.  5-21  that  the  lowest-order  mode,  TEm,  is  excited  at  about  3.38  GHz  which  is 

extremely  close  to  the  resonant  frequency  of  this  mode  (3.40  GHz),  as  shown  in  Table 
5.6.  Also,  the  thin-wire  probe  resonance  is  seen  at  about  2.43  GHz  ,  which  is  close  to  the 
theoretically  [5]  calculated  value  of  2.67  GHz  .  Examining  the  figure,  the  resonances  of 
the  different  slots  are  seen  near  4.48  GHz,  which  is  the  calculated  value.  The  resonance 
of  slot  1  is  seen  at  4.22  GHz  ,  the  resonance  of  slot  2  is  seen  at  4.37  GHz  ,  the  resonance 
associated  with  slot  3  exists  at  about  4.48  GHz  and  the  resonance  of  slot  4  occurs  at 
4.59  GHz .  The  input  impedance  at  the  base  of  the  probe  is  shown  in  Fig.  5-22.  The 
resonances  associated  with  the  cavity,  slot  and  the  probe  are  very  clearly  defined  at  the 
same  frequencies  as  seen  for  the  slot  impedance.  The  EFS  (Fig.  5-23)  is  computed  at  the 
center  of  each  cavity  (z  =  c/ 2,  z  =  3c/2,  z  =  5c/2  and  z  =  7c/2)  using  (5.1).  The 
equation  (5.6)  is  used  to  compute  Ey  at  z  =  7c/ 2 .  The  plot  shows  that  EFS  decreases  in 

the  neighborhood  of  the  resonances  as  expected.  The  EFS  decreases  near  the  vicinity  of 
2.41  GHz  ,  3.40  GHz,  and  4.31  GHz  corresponding  to  the  calculated  resonances. 

The  real  and  imaginary  parts  of  the  input  impedance  at  the  center  of  slot  1  and 
base  of  probe  2  for  Structure  H  in  the  frequency  range  between  6  GHz  and  12  GHz  are 
computed  by  following  the  same  driving  mechanism  as  explained  in  the  previous  section. 
Examining  Fig.  5-24,  the  input  impedance  at  the  center  of  slot  1  is  plotted  and  the  TEm 


mode  is  excited  at  around  4.20  GHz ,  which  is  close  to  the  resonant  frequency  of  this 
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mode  (3.86  GHz).  One  can  also  see  a  slight  change  in  the  input  impedance  at  this 
frequency.  It  can  also  be  seen  from  the  figure  that  the  resonance  associated  with  the  thin- 
wire  probe  occurs  at  about  1.84  GHz,  which  is  close  to  the  theoretical  value  of 
2.08  GHz .  The  resonance  associated  with  slot  1  occurs  at  about  3.02  GHz  which  is 
extremely  close  to  the  calculated  resonance  of  3.0  GHz  .  The  resonance  of  the  second  slot 
is  seen  at  3.07  GHz  while  the  resonances  associated  with  slot  3  and  slot  4  are  not 
properly  excited  in  this  structure.  The  input  impedance  at  the  base  of  the  probe  is  shown 
in  Fig.  5-25  for  a  1.0  V  excitation  at  its  base.  The  resonances  associated  with  the  cavity, 
slot,  and  the  probe  are  very  clearly  defined  at  the  same  frequencies  as  seen  for  the  slot 
impedance.  The  EFS  (Fig.  5-26)  is  computed  at  the  center  of  each  cavity  (z  =  c/ 2, 
z  =  3c/2,  z  =  5c/2  and  z  =  7c/2)  using  (5.1)  for  Structure  H.  Again,  equation  (5.6)  is 
used  to  compute  E  at  z  =  7c/  2 .  The  plot  shows  that  EFS  decreases  in  the  neighborhood 


of  the  resonances,  as  expected. 
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(ju,s) 


(b) 


Fig.  5-19.  Three  cascaded  rectangular  cavities  interconnected  by  narrow  slots  and 
with  thin-wire  probes  in  them,  (a)  perspective  view,  (b)  side  view. 
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Table  5.5 

Structure  Parameters 


Parameter 

Structure  G  (mm) 

Structure  H  (mm) 

Cavity  width,  2 a 

60.0 

53.0 

Cavity  height,  2b 

34.0 

40.0 

Cavity  depth,  c 

65.0 

58.0 

Slot  lengths,  2L,,  2 L2,  2L,,  2L4 

35.0 

50.0 

Slot  width,  2 w 

0.5 

0.5 

Probe  lengths,  /?, ,  h2,  k2 

28.0 

36.0 

Probe  radius,  r 

0.25 

0.25 

Probe  position,  clv  cl2,  d3 

4.0 

20.0 

Table  5.6 

Resonant  Frequencies 


Structure  G 

Structure  H 

TEm 

3.40 

3.86 

te102 

5.25 

5.87 

TEm 

7.36 

7.36 

EE  2oi 

5.5 

5.5 

EE\ 04 

9.56 

9.56 

Probe  half-wave  resonance 

2.67 

2.08 

Slot  half-wave  resonance 

4.28 

3.00 

**  All  Frequencies  are  in  GF 

[z 
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Structure  Parameters 

2Lj  =  2 L,  =  2L3  =  2 L4  =  35.0  mm 

2 w  =  0.5  mm 

1\  =  h2  =  hj  =  24.0  mm 

C  =  65.0  mm 

2a  =  60.0  mm,  2b  =  34.0  mm 
cZ,  =  cl2  =  c/3  =  c/2 


’f  T 


(a)  (b) 


Fig.  5-20.  Magnetic  currents  in  the  slots  and  electric  currents  on  the  probes  due  to  a  1.0 
V  excitation  on  probe  2  (a)  Magnetic  currents  in  slot  1  and  slot  4.  (b) 
Magnetic  currents  in  slot  2  and  slot  3.  (c)  Electric  currents  on  probe  1  and 
probe  3. 
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Fig.  5-21.  Input  impedance  at  the  center  of  slot  1.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Structure  G  Parameters 

2Lj  =  2 L0  =  2 L3  =  2L4  =  35.0  mm 

2w  =  0.5  mm 

=  h2  =  /i,  =  28.0  mm 
c  =  65.0  mm 
dt  =  d1  =  r/3  =  4.0  mm 


Fig.  5-22.  Input  impedance  at  the  base  of  probe  2.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Fig.  5-23.  EFS  for  Structure  G 
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Structure  H  Parameters 

2L,  =  2 L,  =  2L3  =  2 L4  =  50.0  mm 

2w=  0.5  mm 

/7(  =  /j,  =  hj  =  36.0  mm 

c  =  58.0  mm 

=  <j?3  =  20.0  mm 


1  1 
. C  . . wU  c ^ c 


PEC  PEC 


(b) 


Fig.  5-24.  Input  impedance  at  the  center  of  slot  1.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 
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Structure  H  Parameters 
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Fig.  5-25.  Input  impedance  at  the  base  of  probe  2.  (a)  Real  part  of  impedance, 
(b)  Imaginary  part  of  impedance. 


110 


<  > 

c/2 


< -  > 

c+c/2 


<  > 

2c  +  c/2 


< —  > 

3c  +  c/2 


PEC 


PEC 


Fig.  5-26.  EFS  for  Structure  H 
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5.  Conclusions 

In  this  study,  an  analysis  of  thin-wire  probes  inside  slotted  rectangular  cavities  has 
been  undertaken.  The  method  of  analysis  is  based  on  solving  the  coupled  integral 
equations  developed  for  the  electric  current  on  the  probe  and  the  equivalent  magnetic 
current  in  the  slot.  Input  impedance  at  the  center  of  the  slot  and  at  the  base  of  the  probe  is 
calculated  for  single,  double,  and  triply  cascaded  rectangular  cavities.  The  EFS  is 
calculated  at  various  positions  in  the  cavity  to  predict  the  effect  of  the  transmission 
path/environment  on  the  entering  transient  signal.  The  theoretical  calculations  of 
resonances  are  compared  to  the  numerical  results  and  are  found  to  be  in  good  agreement. 
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